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ABSTRACT 

We evaluate pseudoscalar, vector, spin 1/2 and spin 3/2 baryon masses 
predicted by lattice QCD with Wilson quarks in the valence (quenched) 
approximation for a range of different values of lattice spacing, lattice 
volume and quark mass. Extrapolating these results to physical quark 
mass, then to zero lattice spacing and infinite volume we obtain values 
for eight mass ratios. We also determine the zero lattice spacing, infinite 
volume limit of an alternate set of five quantities found without extrapola- 
tion in quark mass. Both sets of predictions differ from the corresponding 
observed values by amounts consistent with the predicted quantities' sta- 
tistical uncertainties. 
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1 INTRODUCTION 



In a recent letter we summarized a lattice QCD calculation of the masses of 
eight low-lying hadrons, extrapolated to physical quark mass, zero lattice spacing, and 
infinite volume, using Wilson quarks in the valence (quenched) approximation. This 
approximation may be viewed as replacing the momentum and frequency dependent 
color dielectric constant arising from quark-antiquark vacuum polarization with its 
zero momentum, zero frequency limit and might be expected to be fairly reliable for 
low-lying baryon and meson masses p[. In the present article we describe the hadron 
mass calculation of Ref. ^ in greater detail. 

Each of the eight predicted mass ratios which we obtain is within 6% of exper- 
iment. The difference between the predicted ratios and experiment is less than 1.6 
times the corresponding statistical uncertainty. It appears to us reasonable to take 
these results as quantitative confirmation of the mass predictions both of QCD and of 
the valence approximation. We believe it is unlikely that the valence approximation 
would agree with experiment for eight different mass ratios yet differ significantly from 
QCD's predictions including the full effect of quark-antiquark vacuum polarization. 

For the four lattices used in our extrapolations to physical mass ratios, we 
calculate hadron masses for a range of quark masses above CSm^, where is the 
strange quark mass. We do not calculate hadron masses directly at smaller quark mass 
because below O.Sm^ our algorithms become unacceptably slow. For quark masses 
ranging from about CSm^ to l.lm^ we find pseudoscalar meson masses squared, vector 
meson masses, and baryon masses to be very close to linear functions of the quark 
mass. The masses of nonstrange hadrons are then found by extrapolating these linear 
fits down to m„, the average of the up and down quark masses. Recent calculations 
based on chiral perturbation theory ||^, ^ show that at sufficiently small quark mass 
a peculiarity of the valence approximation might lead to significant deviations from 
our linear extrapolations. In Ref. p evidence is discussed which suggests that this 
potential difficulty occurs mainly at very small quark mass and probably does not 
have an effect on our extrapolations from 0.3ms down to m„. In the present article we 
will also show that hadron masses in the real world are quite close to linear functions 
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of valence quark masses over the range of our extrapolations. The linearity of observed 
hadron masses as functions of valence quark masses is closely related to the success 
of the Gell-Mann-Okubo mass formula. 

An alternate interpretation of our mass calculations can be made, however, 
which does not depend on extrapolation in quark mass. The linearity of real world 
hadron masses as a function of quark mass implies that all masses in each hadron 
multiplet are determined by the first two coefficients of a Taylor series expansion 
around any quark mass between m„ and m^. Our data shows that for the valence 
approximation this linearity occurs at least between O.Srris and Llm^. Our results 
can then be cast as predictions of the first two Taylor coefficients for each hadron 
multiplet at some conveniently chosen quark mass between O.Srris and Llm^. Four 
of our five predicted values for these coefficients at the point (m„ + ms)/2 differ 
from experiment by less than 1.6 times the corresponding statistically uncertainty. 
The fifth prediction differs by 2.0 times its statistical uncertainty. The predicted 
constant terms are within 6.5% of experiment with statistical uncertainties of up to 
3.3%. The predicted coefficients of the linear terms in these expansions lie within 
22% of experiment with statistical uncertainties of up to 22%. The coefficients of the 
linear terms are obtained from small differences between the masses of hadrons with 
different quark compositions and as a result have relatively large statistical errors. 

Following Refs. P, ^, we also determine the coupling constant from the 
lattice coupling constant giat- From g^^, by the two-loop Callan-Symanzik equation, 
we determine A^|^a in units of lattice spacing a. The values we find for the ratio 
(A^^a) / [rripa) are constant within statistical errors. Thus rripa depends on a accord- 
ing to asymptotic scaling to within statistical errors. This result tends to support the 
reliability of our extrapolation of masses to the continuum limit. The infinite volume, 
continuum limit of A^^/rUp combined with the value of rup in physical units permits 
a calculation of A^^ in physical units. The result we obtain agrees, to within 4% 
statistical errors, with a value of the infinite volume, continuum limit A^^^ obtained 
from a charmonium mass splitting 0. 

An evaluation of meson decay constants using the data set from which masses 
are extracted here is discussed in detail in Ref. |p. 
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The calculations described here were done on the GFll parallel computer at 
IBM Research ||^ and took approximately one year to complete. The machine was 
used in configurations ranging from 384 to 480 processors, with sustained speeds 
ranging from 5 Gfiops to 7 Gfiops. With the present set of improved algorithms and 
480 processors, these calculations could be repeated in about four months. 

2 COULOMB GAUGE HADRON OPERATORS 

For all but one choice of lattice size and (3, we construct hadron propagators 
with nonlocal source and sink operators. The operators which we use are specified in 
lattice Coulomb gauge. A transformation to lattice Coulomb gauge is defined to give 
a local maximum of the sum over all sites and space direction links 

J2 ReTr[U,{x)]. (2.1) 

x,i=l,2,3 

A transformation which produces a local maximum of this sum is found by a method 
qualitatively similar to the Cabbibo-Marinari-Okawa Monte Carlo algorithm. The 
lattice is swept repeatedly, and at each site the target function is maximized first by 
a gauge transformation in the SU{2) subgroup of SU{?>) acting only on gauge index 
values 1 and 2, then by a gauge transformation in the SU{2) subgroup acting only on 
index values 2 and 3, then by a gauge transformation in the subgroup acting only on 
index values 1 and 3. Maximizing the target function over SU (2) subgroups is easier 
to program than a direct maximization over all of S'f/(3). On the other hand, it is 
not clear that maximizing each site over SU{'i) would significantly accelerate the full 
transformation to Coulomb gauge. A local maximum is reached when at each site 
the quantity R{x) vanishes where 

Q{x) = Y.[U,ix)-Ujix~l)], 

i 

R{x) = Q{x) -Q^{x) -hmTr[Q{x)]. (2.2) 

The vector z is a unit lattice vector in the positive i direction. We stop the iteration 
process when the sum over the lattice of the quantity Tr[R\x)R{x)] becames smaller 
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than a convergence parameter c. 

In Coulomb gauge, a smeared quark field 4>r{x, t) is then constructed from the 
local quark field ip{x,t) by 

y 

Gr{z) = (v^r)-3ea;p(-^). (2.3) 

Spin, flavor and color indices have been suppressed in Eq. (p.3|). The field (j)j.{x,t) is 
defined by a corresponding smearing of ipix, t). We take the smeared fields (f)o{x) and 
0q(x) to be ipix) and ^^(x), respectively. 

Smeared hadron fields can then be formed from local products of smeared 
quark and antiquark fields. The fields for a charged pion and a changed rho are 

pt\x) = (2.4) 

For a proton, antiproton, or A++ with with z-component of spin given by s we 
have the nonrelativistic operators 

p:{x) 

Kix) 

a:(x) 

where a, b and c are color indices, i, j, and k are spinor indices, eabc is the alternating 
index, and repeated indices are summed over. For gamma matrices which are the same 
as the Bjorken and Drell convention, but with space direction matrices multiplied by 
i, the nonzero components of the baryon spin wave functions with maximum spin in 
the z-direction are 

pPl/2 _ _pPl/2 _ 
■■-112 ~ 121 ~ 



= 0L(a;)0feir(a;)0ir(a;)eafecrgfc 

= c.(^)0;:,.(^)C.(^)eafecrg^ (2.5) 

= 0air(a;)0L(a;)0L(a;)eafecrgfc 
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pPl/2 _ pPl/2 _ 
J- 443 — 434 — -L) 

rfif = 1, (2.6) 

pA3/2 _ , 
i — i. 

Corresponding smeared fields can be defined for other nucleon and delta charge and 
spin sates, and for other pseudoscalar mesons, vector mesons, and spin 1/2 and 3/2 
baryons. 

Hadron field operators obtained from local products of gaussian smeared fields, 
when fourier transformed, are equivalent to operators formed from local fields with 
nonlocal gaussian relative wave functions. For a pion at rest, for example, we have 

Y,^ns,t)^YlGs{^^iy,th'^|^^y + z,t), (2.7) 

X yz 

where s is v^r. 

Prom Coulomb gauge smeared fields for each hadron we constructed zero mo- 
mentum, spin summed hadron propagators for a collection of source and sink sizes. 
For a hadron /i, the propagator is 

xs 

The invariance of lattice QCD under parity and charge conjugation transformations 
gives a variety of equalities among various pairs of propagators. Propagators which 
should be equal can be added together to decrease statistical fiuctuations. For a 
lattice with time direction period T, our final zero-momentum propagators are then 

C^/(i) = c^j,/(i) -|- c^^i{T — i), 

cut) = <At) + cUT-t), 

C^,,{t) = C,{t) + c^,,{t), (2.9) 
with corresponding definitions for other mesons and baryons. 
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For a lattice with a sufficiently large space direction periodicity S, statistical 
fluctuations in hadron propagators can, in principle, be further decreased by intro- 
ducing several gaussian sources, each multiplied by a random cube root of 1 to cancel 
cross terms between the propagation of different sources for both baryon and meson 
propagators. A closely related idea was ffist proposed in Ref. |T^. If 5* is large and 
even, for example, a quark field with eight sources is 

<f>'r{t) = Y.^yMy,t), (2.10) 



for (j)r of Eq. (|2.3|) . Each component of the eight y included in the sum in Eq. ( p. 10 ) is 
either or S/2, and the for each different y and each different gauge configuration 
used to find quark propagators are independent random cube roots of 1. Using (p'^ 
of Eq. ( p.lO| ) to define hadron source operators but with (pr of Eq. ( p.3| ) still used to 



define sink operators, a new set of hadron propagators C'J^,{t), h = n,p,P,A, can 
be constructed similar to those in Eq. (|2.9| ). If S is sufficiently large, however, the 
C'Jl.,{t) found from any ensemble of gauge fields will behave as though they have been 
averaged over eight times as many independent gauge configurations as corresponding 



3 PROPAGATORS 

Table |I] lists the lattice sizes and parameter values for which quark and hadron 
propagators were evaluated. Up and down quark masses were taken to be equal in all 
propagators. We therefore obtained degenerate masses for each isospin multiplet of 
hadrons. This approximation leads to almost no loss in useful results since since the 
observed values of mass splittings in isospin multiplets in the real world are generally 
somewhat smaller than our statistical errors. We also did no calculations directly 
including strange quarks. Masses for hadrons including strange quarks were found 
by an extrapolation in quark mass to be discussed in Section ^ Direct calculations 
including strange quarks would have required additional programming but would have 
increased our required computer time by only a small fraction. 
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We chose periodic boundary conditions in all directions for both gauge fields 
and quark fields. Gauge configurations were generated using a version of the Cabbibo- 
Marinari-Okawa [0 algorithm adapted for parallel computers by Ding |T^. The 
number of sweeps skipped between configurations and the total count of configurations 
for each set of lattice parameters are given in the last two columns of the table. 
Values of the correlation between hadron propagators on successive pairs of gauge 
configurations we found to be statistically consistent with zero. Thus the number 
of sweeps between configurations in each case was sufficient to produce statistically 
independent hadron propagators on all gauge configurations used in the calculation 
of propagators. 

For the 8^ x 32 lattice at (3 of 5.7, we used point sources and sinks in the 
quark propagators. For all other lattices we used gaussian sources with r of 2, and 
four different gaussian sinks, with r of 0, 1, 2 and 3. For all lattices except 24'^ x 32 
at /5 of 5.7, we found propagators for the source fields of Eq. ( p.3|) including a single 
source. For 24'^ x 32 at /3 of 5.7 we found propagators for the source fields of Eq. ( 2.10 ) 
including eight sources. 

For all the lattices listed in Table |l], the gauge transformation convergence 
parameter c defined following Eq. ( p.2|) was set to 10~^. The average number of 
sweeps required for convergence to this accuracy ranged from 1525 on the lattice 
24^ X 36 at /5 of 5.93, to 2270 on the lattice 16^ x 32 at /3 of 5.7. 

Quark propagators were constructed using the conjugate gradient algorithm 
for the 8^ x 32 lattice at ^ of 5.7, using red-black preconditioned conjugate gradient for 
the other lattices at (3 of 5.7 and 5.93, and using a red-black preconditioned minimum 
residual algorithm at (3 of 6.17 ||13|. At the largest hopping constant values at (3 of 
5.7 and 5.93, preconditioning the conjugate gradient algorithm increased its speed 
by a factor of 3, and at the largest hopping constant values at (3 of 6.17, the change 
from conjugate gradient to the minimum residual algorithm yielded an additional 
factor of 2 in speed. Table gives the average number of lattice sweeps needed to 
find quark propagators for all of the lattice sizes and parameter values which we use 
to find infinite volume continuum limit masses. The number of sweeps in all cases 
was chosen large enough to insure that effective pion, rho, nucleon and delta masses 
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evaluated within the time interval used for final fits are within 0.2% of their values 
obtained on propagators run to machine precision. The number of sweeps required 
to obtain this precision was found by calculating propagators with several different 
convergence criteria on small ensembles of configurations. 

At P of 5.7 we also found propagators for the lattice 16'^ x 32 using a source 
field including eight sources for an ensemble of 81 configurations. For the lattice 
24^ X 32 additional propagators were found using a source field including only a 
single source for an ensemble of 18 configurations. Figures |l| - § show propagator 
statistical dispersions divided by propagators for various combinations of sources and 
lattice sizes. For the lattice 16^ x 32, Figures |1| and ^ compare the dispersion in the 
pion propagators found from a single source and found from eight sources, in both 
cases using 81 gauge field configurations. Figures |^ and § show this comparison for the 
proton propagator. Figures - § compare single source and eight source propagators 
for the lattice 24^ x 32 using 18 gauge field configurations. The range parameter r of 
Eq. (|2.3| ) has the value 2 for all sources and sinks in these figures. The vertical dashed 
lines in each figure mark the fitting intervals, to be discussed later, which we found 
to be optimal for determining the corresponding hadron mass. Figures |l| - § show 
that except for the pion at k of 0.1650, the smallest errors within the fitting intervals 
for propagators on the lattice 16^ x 32 are obtained with a single source. For the 
lattice 24^ x 32, on the other hand, eight sources give significantly smaller statistical 
errors within the fitting interval at k of 0.1650 and at least somewhat smaller errors 
at k of 0.1675. Similar results to those shown in Figures |1| - § for pion and nucleon 
propagators are also obtained for rho and delta baryon propagators. 

Overall it appears that single sources are the most efficient choice for hadron 
propagators for the lattice 16'^ x 32 at /? of 5.7 and therefore also for 24'^ x 36 at /3 
of 5.93 and 30 x 32^ x 40 at 6.17, which have nearly the same volume in physical 
units as 16^ x 32 at /5 of 5.7. For the larger physical volume of 24^ x 32 at /3 of 5.7, 
however, it appears that eight sources are a more efficient choice. 
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4 HADRON MASSES 



Hadron masses were determined by fitting hadron propagators to their asymp- 
totic forms, for large values of t and the lattice time period T. For any meson M and 
any baryon B, 

C^l,{t) Z^^,{exp{-mMt) +exp[-mM{T -t)]}, (4.1) 
C^,,{t) Zf^,exp{-mBt). (4.2) 

Here r and r' are the smearing parameters for propagator sink and source, respec- 
tively. In the asymptotic form for meson propagators we include a contribution from 
backward propagation across the lattice's periodic boundary but omit this term in 
baryon propagators. Backward propagating baryons arising from the operators of 
Eqs. ( p.5| ) are opposite parity excitations with a larger mass than the forward propa- 
gating ground state. For large enough T, with t < T, the omitted term is negligible. 
Including the backward term, on the other hand, would require fitting four parameters 
to baryon propagators rather than two. 

To determine the range of time separations to fit to the asymptotic forms of 



Eqs. and [4.2| , we evaluated for each hadron effective masses m(t) defined to 

be the result of fitting the corresponding propagators to Eqs. ( [4.1D and [4.2| at the 
pair of successive time values t and t + 1. The largest interval at large t showing an 
approximate plateau in an effective mass was chosen as the initial trial fitting range 
for the corresponding propagator. In all cases the initial trial fitting range included 
more than four values of t. An automatic fitting program was then used to choose the 
final fitting range within the initial trial range. For each possible interval including at 
least four values of t within the trial fitting range, the program chose the parameters 
in Eqs. ( [4.1|) and |4.2| which minimize the full correlated of the fit to the data. 



The interval giving the smallest value of P^r fitted degree of freedom and the 
corresponding mass were then chosen as the final fitting range and mass prediction. 
Statistical uncertainties of parameters obtained from fits and of any function of 



these parameters were determined by the bootstrap method |]14| . From each ensemble 



of N gauge configurations, 100 bootstrap ensembles were generated. Each bootstrap 
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ensemble consists of a set of N gauge configurations randomly selected from the 
underlying N member ensemble allowing repeats. For each bootstrap ensemble the 
entire fit was repeated, including a possibly new choice of the final fitting interval. 
The collection of 100 bootstrap ensembles thus yields a collection of 100 values of any 
fitted parameter or any function of any fitted parameter. The statistical uncertainty 
of any parameter is taken to be half the difference between a value which is higher 
than all but 15.9% of the bootstrap values and a value which is lower than all but 
15.9% of the bootstrap values. In the limit of large N, the collection of bootstrap 
values of a parameter p approaches a gaussian distribution and the definition we use 
for statistical uncertainty approaches the dispersion < p'^ > — < p >^. 

In the absence of some independent method for determing the predictions of 
QCD, it appears inevitable that the choice of t interval on which to fit data to a large 
t asymptotic form must be made by some procedure which depends on the Monte 
Carlo data itself. Thus the statistical uncertainties in the data lead to a corresponding 
uncertainty in the choice of fitting interval which, in turn, could lead to an additional 
uncertainty in the fitted result. Another advantage of our procedure for choosing the 
fitting interval combined with bootstrap evaluation of statistical uncertainties is that 
the values we obtain for statistical uncertainties include the uncertainty arising from 
the choice of fitting interval. A comparison of the error bars found for our final fits 
with the error bars found using the same fitting range held fixed across the bootstrap 
ensemble shows that typically about 10% of the final statistical uncertainty comes 
from fiuctuations over the bootstrap ensemble of the fitting range itself. 

Hadron masses were calculated from all propagators discussed in Section |^. 
Figures ^ - 32 show effective hadron masses, fitted hadron masses and final fitting 
ranges for pion, rho, nucleon and delta propagators, for the lattices 16'^ x 32, 24^ x 36 
and 30 x 32^ x 40 which will be used in Section § to obtain continuum limits. Data 
is shown for the lightest quark mass used on each lattice both for sinks with r of 2 
and of 4. Final fitted mass values, fitting ranges and P^r degree of freedom are 



listed in Tables ^ - For all but the lattice 8 x 32, and the lattice 16 x 32 at small 
k, these tables list results for r ranging from to 4 along with masses obtained by 
fitting simultaneously data for sinks with r of 0, 1 and 2. 
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The masses determined simultaneously from sinks 0, 1 and 2 we chose as 
our overall best values. These numbers generally showed somewhat smaller error 
bars than single sink fits. They also provided an unbiased resolution of the small 
disagreements between fits to single sinks of different sizes. We tried simultaneous 
fits to more than three sinks but found that excessive computer time was required 
and the minimal was hard to find reliably. For the pion we found no noticeably 
improvement in errors with multiple sink fits and arbitrarily chose the masses obtained 
from sinks of size as our best values. 

In nearly all cases, the effective masses in the figures show plateaus extending 
over four or more values of t. Since each effective mass is found from a propagator 
at a pair of values of t, it follows that in nearly all cases the propagators shown fall 
off with a consistent mass over five or more values of t. Propagators with r of 4 
show wider plateaus beginning at smaller t than propagators with r of 2. Thus, as 
expected, hadron operators with r of 4 couple more strongly to ground state hadrons 
and less strongly to excited states than do operators with r of 2. The consistency, 
within statistical errors, of mass values shown in Tables |^ - 20 for varying choices of 
r shows that errors in ground state masses which might arise from contamination by 
excited states are within the statistical errors in each mass. 



For larger values of r, the propagators in Figures p| - p2| and the masses in 



Tables - ^ show larger statistical errors. These larger errors arise because hadron 
sink operators with largyer r couple significantly to a larger collection of gauge link 
operators, those between the postions of the sink quark and antiquark fields, and thus 
are more sensitive to fluctuations in the gauge field. 

In subsequent sections of this paper, extrapolations to obtain physical pre- 
dictions for hadron masses will be presented both for mass values obtained from 
simultaneous fits to sinks 0, 1 and 2, and for masses solely from fits to sink size 4. 
The consistency of these two sets of predicitions, within statistical errors, will serve 
as further evidence that our results are not biased by the contamination of ground 
state masses with masses from excited states. 
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5 EXTRAPOLATION TO SMALL QUARK MASS 



As the hopping constant k is made larger, the amount of computer time re- 
quired to evaluate a quark propagator grows, as does the size of the ensemble of gauge 
field configurations required to find hadron propagators to within a fixed statistical 
error. As a consequence of these effects, particularly the second, we were unable to 
evaluate hadron masses at large enough k to give the physical value of 171^,/ nip. Thus 
to determine the masses of hadrons composed only of up and down quarks we ex- 
trapolated to the physical value of k data found at smaller k on the lattice 16'^ x 32, 
24^ X 32, 243 X 36 and 30 x 32^ x 40. For the lattice 8^ x 32 we were not able to run 

ge enough to make this extrapolation reliably. 

The quark mass in lattice units m^a can be defined as 

rriqa = — — , (5.1) 

where k^ is the critical hopping constant at which becomes zero. A naive applica- 
tion of the hypothesis that the pion is the pseudo-Goldstone boson of spontaneously 
broken chiral symmetry implies that should be linear in rriq if rriq is made small 
enough. Eq. ( ^.1|) then yields linearity in l/{2k) near l/(2A;c). On the lattices 16^ x 32, 
24^ X 32, 24^ X 36 and 30 x 32^ x 40, we found ml to be close to a linear function of 
1 / {2k) at all the values of k in Table 0, and statistically consistent with exact linearity 
at the three largest values of k. The best fit at the three largest values of k in each 
case was found by minimizing the fit's full correlated x^- Values of k^ were taken from 
these fits. Eq. (|5.1| ) then gives a translation from k to rrig. For rrip, m^, and rriA, 
a simple application of perturbation theory suggests linearity in ruq at small enough 
TTiq. These quantities we found to be close to linear in rriq at all k and statistically 
consistent with exact linearity at the three largest k. Figures ^ - ^ show these fits 
and extrapolations for masses from sinks 0, 1 and 2, and for masses from sink 4 for 
the lattices 16^ x 32, 24^ x 32, 24^ x 36 and 30 x 32^ x 40. Hadron masses in Figures || 
- 38| are shown in units of nip extrapolated to physical quark mass, and the quark 
mass is shown in units of the strange quark mass nig. How we determine will be 
discussed in Section ^. 
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On each lattice, the fits of and rrip to hnear functions of m^a we used to 
determine the value of m^a which gives the physical value of m^/mp. This choice of 
rUqa we define to be the "normal quark" mass m„a. The fits, for each lattice, of rriNCi 
and rriACi to linear functions of m^a we then evaluated at m^a to find mAra and mAfl 
at physical quark mass. These masses taken at m„a in all cases differ by less than one 
standard deviation from their values extrapolated to zero ruq. Numbers for m„a and 
''Tip{mn)a are given in Table 21. Tables give a variety of ratios of extrapolated 

hadron mass values both for data obtained from simultaneous fits to sinks of size 0, 
1 and 2 and for data from fits to sinks of size 4. For each ratio in each table, the last 
column gives the value of per degree of freedom of the linear fit used to extrapolate 
the ratio's numerator to the correct quark mass. 

Following Refs. [0, P] we have also calculated, for each lattice, values of g^^- 
From these and the two-loop beta function for the Callan-Symanzik equation, we 

of mp{mn)a. 



have found for each lattice A^-j-a. Tables 22 - 25 give these results measured in units 

i W O ' ' ' ' 



6 STRANGE HADRON MASSES 

The linear relations between hadron masses and quark and antiquark masses 
which we have found to be approached at small enough quark and antiquark mass 
for hadrons composed of a single mass of quark and antiquark can be summarized 

m^{mqY = a^rriq, (6.1) 
mp{mq) = apTUq + Cp, (6.2) 
rriBimq) = asmq + CB- (6.3) 

Here rriq is a variable mass put in place of the normal quark and antiquark mass, and 
B is either the nucleon or the delta baryon. These equations strongly suggest that for 
strange hadrons composed of quarks and antiquarks with different masses we should 
have in addition 

mK{mqi,mq2y = axrUqi + bKmq2, (6.4) 
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mK*{rnqi,mq2) = aK*mgi + bK*mq2 + Ck*, (6.5) 
rriBimgi, mq2) = aBirigi + bBmq2 + Cb, (6.6) 

where m^i and are variable masses put in place of the normal and strange quark 
masses, respectively, and B is any strange baryon. 

Assuming Eqs. ( |6.4| ) - ( |6.6| ), a set of equations can be derived among the coef- 
ficients of related hadrons. These equations permit the masses of hadrons containing 
both strange quarks and normal quarks to be obtained from the masses of related 
hadrons composed of a single flavor of quark with quark mass some weighted average 
of normal and strange quark masses. The equations among masses found in this way 
are closely connected to the Gell-Mann-Okubo mass formulas. 

For the pseudo-scalar and vector mesons, charge conjugation invariance implies 

ax = bx, 

ax* = bx*- (6.7) 

By flavor SU(3) symmetry, pion and kaon masses become degenerate if rriqi and mq2 
are equal, and the rho and k-star masses become equal if rriqi and 171^2 are equal. 
Thus 

ax = CL-K, 
ax* = ctp, 

CX ^TT 5 

cx* = Cp. (6.8) 

We obtain for strange meson masses at physical values of the normal and 
strange quark masses 

mx{mn,ms) = m^[(m„ + mJ/2], (6.9) 

mx*{mn,ms) = mp[{mn + ms)/2]. (6.10) 

In addition, in the valence approximation, the phi vector meson does not mix with 
the omega and is composed purely of a strange quark and antiquark, giving for the 
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physical phi mass 



(6.11) 



For baryons, relations similar to Eqs. ( |6.9D - ( |6.11| ) can be derived by combining 
Eq. (|6.3| ) and (|6.6| ) with the asymptotic form Eq. ( [4.2| ) and definitions Eqs. ( |2.5| ) - 
(|2.6|) . Differentiating Eq. ( [4.2|) with respect to a quark mass and taking the asymptotic 
behavior for large t gives 



d 



dm, 



-C^{t) 



qi 



d 



dm, 



-rriB 



(6.12) 



Evaluating the left side of Eq. ( |6.12| ) for m^i equal to mq2 and using Eqs. (|2.5|) - (p.6|), 
a variety of linear relations can be obtained among the derivatives of baryon masses 
with respect to quark masses. 

For the spin 1/2 baryon multiplet we find 



a-N = bj^ + bs, (6.13) 

bs = lbj: + h^, (6.14) 

and the coefficients cb are the same for all members of the multiplet. For the spin 
3/2 baryon multiplet we find 



bn 
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OA, 



(6.15) 

(6.16) 
(6.17) 



and the coefficients cb are, again, the same for all members of the multiplet. 
From Eqs. ( |6.13|) - (|6.17|) we obtain 



m-sinin, ms) + ms(m„, m^) - mN{mn) 

mj:,{mn,ms) 
ms* {mn,ms) 
mnims) 



mA[(2m„ + ms)/3], 
mA[{mn + 2mJ/3], 
mA(mJ. 



(6.18) 
(6.19) 
(6.20) 
(6.21) 
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Eqs. ( |6.9| ) - ( |6.11| ) and ( |6.18| ) - ( |6.21| ) permit a variety of strange hadron masses 



to be obtained from the fits in Section ^ of hadron masses to Eqs. ( |6.1| ) - (|6.3D. For 
each lattice, we determine the strange quark mass by tuning rris in Eq. ( |6.9|) to 
produce the physical value of ttik/^p- For the values of m„ and irtp needed in these 
calculations we take the results of Section |^. Values of found in this way are 



listed in Table 21. Strange hadron masses determined from m„ and are given 



in Tables 22-05^ both for simultaneous fits to sinks of size 0, 1 and 2 and for fits 



to sink size 4. The hadron masses in these tables are all measured in units of the 



extrapolated nip given in Table |2l|. For each mass ratio in each table, the last column 
again gives per degree of freedom for the linear fits used to determine the ratios 
numerator. 

7 TEST OF EXTRAPOLATION TO SMALL QUARK MASS 



A test of the linear relations Eqs. (|6.1|) - (|6.6|) and of our determination in 
Section ^ of the masses of light hadrons by extrapolation can be made using observed 
values of hadron masses. 

Eqs. (16.101) and (|6.11| ) combined with physical values of the k-star and phi 



masses give values for the mass of a rho made of heavy quarks. Eqs. ( |6.19D - (|6.21|) 
can be used similarly to find the mass of a delta baryon composed of heavy quarks. 
For the nucleon, define m^A to be 



Then we have 



ms(m„,ms) = msA("^n, 2ms), (7.2) 
rriNimn) = mY.K{mn,mn). (7.3) 



The first of these equations follows from Eq. ( |6.14| ) , and the second holds because the 



masses of the octet of spin 1/2 baryons becomes degenerate if irts and are equal. 

From mp, m/\ and m^A for heavy quarks we can now attempt to recover the 
physical values mp{mn), iTi/^{mn) and m7v(m„) by linear extrapolation as done in 
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Section ^ Figure ^ shows linear extrapolations in niq down to the value m„ of 
mpiniq), mY;A{mn, rUq) and m/^{mq). The scale for is shown in units of the strange 
quark mass and the hadron mass scale is shown in units of the physical rho mass 
^p(^n)- The hadron mass at equal to m„ compared with the msA extrapolation, 
following Eq. ([7.3|), is the physical nucleon mass mj^{mn)- The extrapolated value for 
^pi^n) is low by 0.53%, for mN{mn) is high by 1.38%, and for m^{mn) is high by 
0.81%. The linear fit in each of these extrapolations uses a range of mq/rris contained 
within the range used in our extrapolations in Section |^. Overall these results support 
the accuracy of the linearity assumed in Section ^ and the extrapolation to find hadron 
masses in Section |. 



8 CONTINUUM LIMIT 

The value of p for each of the three lattices 16^ x 32, 24^ x 36 and 30 x 32^ x 40 
was chosen so that the physical volume in each case is nearly the same. For lattice 
period L, the quantity rUpL is respectively, 9.08 ± 0.13, 9.24 ± 0.19 and, averaged 
over three directions, 8.67 ± 0.12. Thus a sequence of corresponding results on these 
three lattices gives each predicition's behavior as the lattice spacing is made smaller 
with the volume held fixed in physical units. 

For sufficiently small values of lattice spacing, the leading lattice spacing de- 
pendence of hadron mass ratios is expected to be linear in a. Figures ^ - ^ show 



mass ratios for the lattices 16^ x 32, 24^ x 36 and 30 x 32^ x 40, listed in Tables ||, |4 



and as a function of lattice spacing measured in units of mp(m„). The lines are 
linear fits to these points. The vertical bars at zero lattice spacing are the statistical 
uncertainties in the linear extrapolations of hadron mass ratios to zero lattice spacing. 
The dots at zero lattice spacing represent the observed physical values of ratios. The 
extrapolated ratios, corresponding observed values and per degree of freedom of 
the linear fits are given in Table |2^. 

The finite volume continuum limits shown in Table ^for simultaneous fits to 
sinks 0, 1 and 2 and for fits to sink 4 are consistent and all lie within 1 standard devi- 
ation of each other. The size of the statistical errors for the simultaneous fits to sinks 
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0, 1 and 2 are smaller than those for sink 4. In general, for an increasing sequence of 
Monte Carlo statistical uncertainties, the uncertainties in these uncertainties increase 
more rapidly. Thus we believe the simultaneous fits to sinks 0, 1 and 2 and their error 
bars are more reliable than the corresponding numbers found from sink 4. With the 
exception of the value of (m^ + — m^)/mp, the predicted numbers lie within 1.7 
standard deviations of the observed results and are statistically consistent with the 
observed results. The prediction for (ms + ttlj: — mN)/mp differs from experiment 
by 4.0 standard deviations. We will show in Section ^, however, that the agreement 
between the prediction for (m= + mE— mAr)/mp and observation becomes much better 
when a correction to obtain infinite volume results is applied. 

Figures ^ and |47| show linear extrapolations of Apr^/mp(m„) to the continuum 
limit, for mp{mn) determined from sinks 0, 1, 2 and from sink 4, respectively. The 
dot at zero lattice spacing in these figures is the value determined from heavy quark 
spectroscopy in Ref. [H]. Numerical values of the continuum limits shown in these 



figures are listed in Table 26 along with the result from Ref. 



In both Figures 56 and J7, the slope of the linear fit is statistically consistent 
with zero. Thus it follows that the rho mass in lattice units mp{mn)a depends on 
as predicted by the Callan-Symanzik equation using the two-loop beta function. 
Figure shows mp{mn)a, for sinks 0, 1, 2, as a function of c^jj^ for the lattices 
16^ X 32, 24^ X 36 and 30 x 32^ x 40. The line in this figure is the prediction of the 
Callan-Synamzik equation with A^^/mp(m„) given by the continuum limit value in 



Table ^ found by linear extrapolation. The data in Figure ^ appears to be quite 
close to the asymptotic scaling curve. This evidence for asymptotic scaling tends to 
support the reliability of the linear extrapolations we have used to find continuum 
mass ratios. 

The fits in Figures [4^ - ^ were done by minimizing the from the full 
correlation matrix among the fitted data. Both the x and y coordinates of each of the 
three fitted points on each line have statistical uncertainties. We therefore evaluated 
among all six pieces of data and chose as fitting parameters the slope and intercept 
of the line along with the x coordinate of each point. The correlation matrices which 
we used were found by the bootstrap method as were the statistical uncertainties of 
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the extrapolated predictions. The correlation matrices used in fits for each bootstrap 
ensemble were taken, for convenience, from the full ensemble and not recalculated on 
each bootstrap ensemble independently. 



9 INFINITE VOLUME LIMIT 

The continuum limits found so far are for a finite volume lattice. We now 
consider the infinite volume limit of these finite volume continuum predictions. 

As a first step, we compare masses found on the lattices 16^ x 32 and 24^ x 32 



at P of 5.70. For the values of k at which we did direct mass calculations. Table ^ 
shows percent changes in masses from 16^ x 32 to 24^ x 32 at /? of 5.70. Data is 
shown for simultaneous fits to sinks 0, 1 and 2 and for fits to sink 4. For masses from 
fits to sinks 0, 1 and 2 there is some indication, with marginal statistical significance, 
of decreases of up to about 5%. For masses from fits to sink 4, there is still weaker 
evidence for decreases in mass ranging up to about 3%. Overall, it appears to us 



the data in Table ^ is best taken as evidence for an upper bound of 5% on the 
change in masses from 16'^ x 32 to 24^ x 32 at /? of 5.70 for the values of k at which 
we calculated hadron masses directly. For hadron masses extrapolated to physical 
quark mass, measured in units of mp(m„). Table provides a bound of about 
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on volume dependence both for mass ratios from simultaneous fits to sinks 0, 1 and 
2 and for masses from fits to sink 4. 

From the change between a hadron's mass evaluated on a lattice 16^ x 32 and 
on a lattice 24^ x 32, an estimate can be made of the change from 24^ x 32 to infinite 
volume if we have sufficient information concerning the form of the dependence of 
hadron masses on lattice volume. A simple non-relativistic potential model implies 
that the error in a particle's mass due to calculation in a finite volume falls at 
large L as Ce~^^^, where R is the radius of the probability density of the hadron's 
wave function and C is either constant or monotonically decreasing. At f3 of 5.70, R 
for each of the hadrons we consider is then less than about 4 lattice units. On the 



other hand, a model |[T^ based, in part, on a rigorous argument |]T6| gives a volume 
dependent error of the form DL^^ with a constant D for the range of L we consider. 
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At still larger L this model yields an error falling exponentially. 

Assuming volume dependent errors of the form Ce~^^^, the changes we have 
found between masses on a lattice 16^ x 32 and on a lattice 24^ x 32 imply changes in 
mass from 24^ x 32 to infinite volume of less than 1%. Assuming volume dependent 
errors of the form DL~^, the changes we have found between 16^ x 32 and 24^ x 32 
imply changes in mass from 24^ x 32 to infinite volume of less than 2%. 

Estimates can now be made of the corrections needed to obtain infinite volume 
continuum mass ratios from our finite volume values. For the ratio of any hadron 
mass to the rho mass nih/mp, define the finite volume correction term 6h{a, L) to be 

Sh{a,L) = — {a, oo) -—{a, L). (9.1) 
nip nip 

The quantity which we would like to determine is 5/^(0, 9/mp), since 9/mp is L for 
the lattices we used to find continuum limit masses. In Section ^ nih/nip for all h we 
considered shows a relative change of less than 20% as a goes from its value 05.7 at 
(3 of 5.7 down to 0. Thus we would expect a corresponding error of less than 20% of 
(5/i(0,9/mp) for the approximation 

5,(0,— )^(5,(a5.7,— ). (9.2) 
nip nip 

On the other hand, from our estimate of the volume dependent error in masses found 
on 24^ X 32, for which L is 13.5/mp, it follows that with an additional of less than 
1% or 2% of nih/nip we have 

r , 9 nih 13.5 nih, 9 

^h{a5.7, — ) ~ — (^5.7, ) (05.7, — ), (9-3) 

nip nip nip nip nip 

Finally, for each of the hadrons we consider, we have already found that the right 



side of Eq. ( |9.3|) is less than 5% of nih/nip. Combining Eqs. ( pTll) - (|9.3| ), we obtain 
the approximation 

ruh,^ \ ^hf^ 9 nih, 13.5 nih 9 

— (0, 00) ^ — (0, — ) + — (05.7, ) (05.7, — ), (9.4) 

rrir, nin nir, nin nir, nir, m„ 
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with two contributions to the error. One contribution is less than 1% or 2% of nih/mp, 
and the other less than 20% of 5% of mh/mp, which is 1% of mh/rrip. Thus overall 
the error in Eq. ( |9.4D should be less than about 2% of rrih/mp. 

Table |2^ shows infinite volume continuum limit mass ratios obtained from 
the finite volume values using Eq. ( |9.4|) . Results again are shown both for masses 
from simultaneous fits to sinks 0, 1 and 2 and from fits to sink 4. The two sets of 
numbers are statistically consistent but the errors for the fits to sinks 0, 1 and 2 are 
smaller. As before, we consider both the predicted values for sinks 0, 1 and 2 and 
the statistical errors on these numbers to be more reliable than the corresponding 
numbers for sink 4. The main effect of the correction to infinite volume is an increase 
in standard deviations by about a factor of 1.5. The shift in central values in all 
cases is small and less than about 1.2 infinite volume standard deviations and 1.9 
finite volume standard deviations. The only significant effect of the infinite volume 
correction on the comparison between predicted numbers and experiment is for the 
value of {m-^ + mY_ — m]^)/'mp. The finite volume prediction for sinks 0, 1 and 2 differs 
from experiment by 4.0 standard deviations while the infinite volume number differs 
from experiment by only 1.6 standard deviations. 

The infinite volume continuum limit predictions for sinks 0, 1 and 2 in Ta- 
ble ^ are statistically consistent with experiment. The predicted values differ from 
experiment by amounts ranging up to 1.6 standard deviations. As a fraction of the 
observed results, the errors range up to 6% with statistical uncertainties ranging up 
to 8%. The infinite volume continuum limit predictions for sink 4 are also statisti- 
cally consistent with experiment but do not agree with experiment as well as do the 
predictions from sinks 0, 1 and 2. The errors for the sink 4 predictions range up to 1.6 
standard deviations, with the exception of a single error of 2.5 standard deviations. 
As a fraction of observed results, the errors go up to 10% with uncertainties up to 
11%. In place of an observed value for A^^^ with which to compare our corresponding 
prediction in Table |2^, we use an infinite volume continuum limit result obtained by 
a lattice QCD calculation combined with the observed value of a charmonium mass 
splitting p. 
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10 PREDICTIONS WITHOUT EXTRAPOLATION TO SMALL QUARK MASS 



Calculations based on chiral perturbation theory p, ^ show that at sufficiently 
small quark mass a peculiarity of the valence approximation might lead to significant 
deviations from the linearity between hadron masses and quark mass which we find 
for quark masses above 0.3ms. In Ref. evidence is discussed which suggests that 
this potential difficulty occurs primarily at very small quark mass and probably does 
not have a significant effect on our extrapolations from O.Srris down to rrin- We now 
consider, however, an alternate interpretation of our results which does not depend on 
the extrapolation of hadron masses beyond the interval within which we have direct 
evidence for linearity. 

The linearity of real world hadron masses as a function of quark mass found in 
Section |^ implies that all masses in each hadron multiplet are determined by the first 
two coefficients of a Taylor series expansion around any quark mass between m„ and 
rris. A convenient expansion point, for example, is rrins defined to be (m„ + ms)/2. We 
now reanalyze our data to obtain predictions for five of the eight significant coefficients 
in Taylor expansions of mT,{mg), mp{mq), mN{mq) and mA{rnq) as functions of rrig — 
irtns- The three coefficients not predicted are, in effect, taken from experiment and 
used to fix the three free parameters of lattice QCD. 

For each of the four lattices 16^ x 32, 24^ x 32, 24^ x 36 and 30 x 32^ x 40, the 
fits in Section]^ of mT,{mq)^, nip^mg), mj^irriq) and m/^{mq) to linear functions oiniqa 
we interpret, using Eq. |5.1| , as fits to linear functions of l/(2fc). We thereby avoid the 
implicit dependence on extrapolation built into the original fits as a consequence of 
the definition of m^a by Eq. (|5.1|) . The critical hopping constant kc entering Eq. ( |5.1| ) 
is found, in effect, by extrapolation. We then determine the hopping constant kns 
corresponding to m„s by requiring mT,{mns) / pij^ns) to agree with the real world 
value of raxIrnK* as expected according to Eqs. ( |6.9|) and ( |6.10| ). For any k we define 
the quark mass difference 

TJiqa - rrinsa = (10.1) 

The hopping constant /c2ns corresponding to 2m„s we find by requiring mT^{2mns) / ^p{fn 
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to be equal to ^J^mxIvriK* as expected according to Eq. (|6.1|). Both kns and k2ns for 
all four lattices lie within the interval for which the fits in Section ^ give direct evi- 
dence, without extrapolation, of linearity of m.„{mgY, mp{mq), mM{mq) and miy{mq) 
in l/(2fc). The determination of kns and k2ns and the determination of hadron masses 
at these points requires only interpolation of our data, not extrapolation. The value 
in lattice units m„sa is found from Eq. (|10.1|) to be 

rrinsa = ^ - (10-2) 

A definition of m^a which does not depend on extrapolation can be found by combing 
Eqs. ( [10. 1| ) and ( |10.2| ). For the present discussion, however, we do not need a definition 
of m^a itself but only the difference m^a — m^^a. 

Our fits m^(mg)^, mp{mq), mNlmq) and 'm/\{mq) to linear functions of l/{2k) 
we now reinterpret, using Eq. ( |10.1| ), as fits to linear functions of m^a — rrinsa. Of the 
eight coefficients entering these fits, two have been chosen in the course of determining 
kns and fc2ns- A third is needed to determine the lattice scale. We are left with five 
predictions of Taylor coefficients at the point at which rriqa — mnsd is 0. For the four 
lattices 16^ x 32, 24^ x 32, 24^ x 36 and 30 x 32^ x 40, these five predictions are shown 
in Tables ^ - 13^. We show also predictions for Aj^/mp{mns)- Numbers are again 
given both for fits to sinks 0, 1 and 2, and for fits to sink 4. The two sets of numbers 
are statistically consistent, but the error bars for the fits to sink 4 are significantly 
larger than those for the fits to sinks 0, 1 and 2. The last column in each table gives 
the degree of freedom of the corresponding linear fit from which each mass or 

mass derivative was determined. 



Table ^ shows the continuum limits of these predictions taken with physical 
volume held fixed, the continuum limits corrected to infinite volume, and the corre- 
sponding observed values. The continuum limit predictions are found as in Section |^ 
by linear extrapolation to zero lattice spacing. We now use mpirnns)^ as the measure 
of lattice spacing, rather than mp{mn)a,. The per degree of freedom for each linear 
fit used to find a continuum limit with physical volume held fixed is shown in the 
last column of Table Figures ^ and ^ show the extrapolations to zero lattice 
spacing with physical volume held fixed. The infinite volume predictions shown in 
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Table ^ are found following Section ^ Again the results for simultaneous fits to sinks 
0, 1 and 2 are statistically consistent with the results for sink 4, but the statistical 
uncertainties for the fits to sink 4 are significantly larger than those for fits to 0, 1 
and 2. We consider the 0, 1, 2 predictions to be more reliable. 

Observed numbers in Table |3^ are determined as discussed in Section 0. The 
infinite volume continuum limits of the predictions for sinks 0, 1 and 2 are statistically 
consistent with the corresponding observed values. Four of the five predicted values 
are within 1.6 standard deviations of experiment and the fifth differs from experiment 
by less than 2.0 standard deviations. The mass predictions are within about 6.5% 
of experiment with statistical uncertainties of up to 3.3% of experiment. The slope 
predictions are within 22% of experiment with statistical uncertainties of up to 22% 
of experiment. The slope predictions are obtained, in effect, from comparatively small 
differences between predicted masses for hadrons composed of quarks with different 
masses and therefore have larger relative statistical errors than the mass predictions. 
As in Table 26, in place of an observed value for A^^ in Table 33, we use an infinite 
volume continuum limit result obtained by a lattice QCD calculation combined with 
the observed value of a charmonium mass splitting [0]. 
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Table 1: Number of configurations analyzed for various lattices sizes and 
Quark masses niq are given in MeV/c^. 
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lattice (3 algorithm k sweeps 



16^ X 32 5.70 conj. grad. 0.1600 500 

0.1650 1000 

0.16625 1500 

0.1675 4000 

24^ X 32 5.70 precond. 0.1600 150 

conj. grad. 0.1650 340 

0.1663 511 

0.1675 998 

24=^ X 36 5.93 precond. 0.1543 195 

conj. grad. 0.1560 240 

0.1573 380 

0.1581 724 

32^ X 30 X 40 6.17 precond. 0.1500 184 

min. resid. 0.1519 350 

0.1526 709 

0.1532 1389 

Table 2: Algorithms and average number of sweeps used to find quark propagators. 
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11 


- 14 


0.2 




0.1500 





1.8416 ±0.0053 


11 


- 14 


0.6 




0.1550 





1.6146 ±0.0076 


12 


- 15 


0.5 




0.1600 





1.3959 ±0.0159 


10 


- 13 


0.3 




0.1650 





1.2476 ± 0.0337 


7- 


- 10 


2.4 



Table 3: Masses from an 8^ x 32 lattice at (3 of 5.70 using 2439 configurations. 



29 



particle 


k 


sink 


ma 


t 




TT 


0.1400 





1.4109 ±0.0057 


10-13 


0.0 




0.1450 





1.2423 ±0.0052 


10-13 


0.0 




0.1500 





1.0691 ±0.0047 


10-13 


0.0 




0.1550 





0.8891 ±0.0040 


9 -12 


0.1 


P 


0.1400 





1.4351 ± 0.0072 


10-13 


0.0 




0.1450 





1.2761 ± 0.0068 


10-13 


0.0 




0.1500 





1.1183 ±0.0064 


9 -13 


0.1 




0.1550 





0.9562 ± 0.0066 


12-15 


0.2 


N 


0.1400 





2.2763 ±0.0235 


8 -11 


1.0 




0.1450 





2.0307 ±0.0207 


8 -11 


0.7 




0.1500 





1.7844 ±0.0205 


8 -11 


0.4 




0.1550 





1.5358 ±0.0227 


8 -11 


0.1 


A 


0.1400 





2.2765 ±0.0252 


11-14 


1.1 




0.1450 





2.0473 ± 0.0244 


8 -11 


0.8 




0.1500 





1.8088 ±0.0227 


8 -11 


0.6 




0.1550 





1.5736 ±0.0256 


8 -11 


0.4 



Table 4: Masses from a 16^ x 32 lattice at /5 of 5.70 using 47 configurations. 



30 



k 


sink 




t 






0.1600 





0.6889 ±0.0013 


11 - 


14 


0.2 




1 


0.6891 ±0.0015 


10 - 


14 


0.2 




2 


0.6904 ±0.0021 


10 - 


13 


0.4 




3 


0.6917 ±0.0028 


10 - 


13 


0.3 




4 


0.6929 ±0.0032 


10 - 


13 


0.2 


0.1650 





0.4584 ±0.0018 


11 - 


14 


1.0 




1 


0.4590 ±0.0019 


10 - 


13 


0.4 




2 


0.4601 ± 0.0021 


lo- 


13 


0.3 




3 


0.4582 ± 0.0030 


ll - 


14 


0.5 




4 


0.4579 ± 0.0034 


11 - 


14 


0.2 


0.16625 





0.3866 ±0.0022 


11 - 


14 


0.5 




1 


0.3869 ±0.0021 


10 - 


13 


0.4 




2 


0.3880 ±0.0025 


10 - 


13 


0.2 




3 


0.3865 ± 0.0029 


11 - 


14 


0.3 




4 


0.3865 ± 0.0033 


11 - 


14 


0.1 


0.1675 





0.2977 ±0.0043 


11 - 


14 


0.1 




1 


0.2979 ± 0.0040 


11 - 


14 


0.1 




2 


0.2979 ± 0.0040 


10 - 


14 


0.1 




3 


0.2981 ± 0.0046 


11 - 


14 


0.2 




4 


0.2991 ± 0.0051 


11 - 


14 


0.1 



Table 5: Values of m^ra from a 16'^ x 32 lattice at jS of 5.70 using 219 configurations. 



31 



k 


sink 


TTlpd 




t 


A 


0.1600 





8044 ± 0029 


9 ■ 


- 13 


0.5 




1 


8061 ± 0030 


10 


- 13 


0.0 




2 


8090 ± 0035 


10 


- 13 


0.0 




3 


8113 ± 0058 


11 


- 14 


0.0 




4 


8133 ± 0070 


11 


- 14 


0.1 


0.1650 


012 


6642 ± 0052 


7 . 


- 10 


0.9 







0.6621 ± 0.0062 


7 ■ 


- 11 


0.8 




1 


6660 ± 0050 


6 


- 9 


0.5 




2 


6627 ± 0054 


6 


- 9 


0.0 




3 


6598 ± 0074 


6 


- 9 


0.2 




4 


6705 ± 0070 


2 


- 6 


0.7 


1 6625 


012 


6321 + 0062 


7 . 


- 10 


4 







0.6306 ±0.0066 


7 ■ 


- 10 


0.3 




1 


0.6302 ±0.0064 


7 ■ 


- 10 


0.2 




2 


0.6311 ±0.0063 


6 


- 9 


0.0 




3 


0.6281 ±0.0080 


6 


- 9 


0.1 




4 


0.6262 ±0.0087 


4 


- 8 


0.4 


0.1675 


012 


0.6068 ±0.0051 


5 


- 8 


0.4 







0.6050 ± 0.0062 


5 


- 8 


0.1 




1 


0.6053 ±0.0061 


4 


- 7 


0.2 




2 


0.6017 ±0.0075 


5 


- 8 


0.3 




3 


0.5940 ±0.0085 


4 


- 7 


0.8 




4 


0.6006 ±0.0077 


2 


- 6 


0.7 



Table 6: Values of m^a from a 16^ x 32 lattice at (3 of 5.70 using 219 configurations. 



32 



k 


sink 


mj\ra 




t 


x' 


0.1600 





1.3043 ±0.0080 


8 ■ 


- 12 


0.6 




1 


1.3015 ±0.0077 


7- 


- 10 


0.2 




2 


1.3004 ±0.0090 


6 ■ 


- 10 


0.1 




3 


1.2905 ±0.0162 


6 


- 9 


0.4 




4 


1.2985 ±0.0141 


2 


- 5 


0.8 


0.1650 


012 


1.0312 ±0.0125 


8 ■ 


- 11 


0.8 







1.0276 ±0.0143 


8 ■ 


- 11 


0.9 




1 


1.0293 ±0.0182 


7 ■ 


- 10 


0.7 




2 


1.0293 ±0.0176 


6 ■ 


- 10 


0.7 




3 


1.0419 ±0.0198 


3 


- 8 


0.7 




4 


1.0322 ±0.0154 


1 


- 4 


0.1 


0.16625 


012 


0.9838 ±0.0103 


5 


- 8 


1.2 







0.9654 ±0.0138 


7- 


- 10 


1.4 




1 


0.9867 ±0.0188 


4 


- 7 


0.3 




2 


0.9857 ±0.0148 


3 


- 6 


0.4 




3 


0.9708 ±0.0124 


3 


- 8 


0.4 




4 


0.9576 ±0.0178 


1 


- 4 


0.1 


0.1675 


012 


0.9052 ± 0.0123 


5 


- 8 


0.5 







0.9100 ±0.0115 


5 


- 8 


0.1 




1 


0.9120 ±0.0110 


5 


- 8 


0.0 




2 


0.9065 ±0.0124 


3 


- 7 


0.3 




3 


0.8906 ±0.0128 


2 


- 5 


0.3 




4 


0.8610 ±0.0154 


2 


- 5 


0.1 



Table 7: Values of rriNa from a 16^ x 32 lattice at (5 of 5.70 using 219 configurations. 



33 





OllllV 


777 A n 




u 


A 


0.1600 





1 3804 ± 0088 


7 . 


- 12 


1.4 




1 


1 371 3 + n 0088 


7 . 

f 


- 10 


4 




2 


1 3687 ± 0173 


6 ■ 


- 10 


3 






1 3Q80 ± 0130 


2 


- 6 


1 




4 


1 3703 + n 014Q 


2 


- 5 


5 




019 






- 1 1 

" X X 


1 







1 1597 ± 0132 


7 . 


- 10 


1.1 




1 


1 1519 ± 0268 


7 . 


- 10 


0.5 




2 


1 1 845 + 0261 


4 


_ 7 


1 

X . \J 






1 1876 + 014Q 


2 


- 6 


6 




4 


1 1511 ± 0211 


1 


. 5 


0.3 


1 fifi25 

yj . ± yjyj^fj 


01 9 


1 1 91 4 + n ni 80 

X . X X t: I VJ . VJ X (J W 


7 . 


- 1 

X w 


1 4 







1.1119 ± 0.0144 


7 ■ 


- 10 


0.4 




1 


1.1054 ±0.0273 


7 ■ 


- 10 


0.2 




2 


1.1520 ±0.0171 


3 


- 6 


0.5 




3 


1.1481 ±0.0183 


1 


- 4 


0.3 




4 


1.1012 ±0.0215 


1 


- 5 


0.3 


0.1675 


012 


1.1349 ±0.0159 


5 


- 8 


2.1 







1.1208 ±0.0141 


5 


- 8 


1.0 




1 


1.1164 ±0.0134 


4 


- 7 


0.2 




2 


1.1067± 0.0141 


4 


- 7 


0.1 




3 


1.0874 ±0.0170 


2 


- 5 


0.1 




4 


1.0526 ±0.0167 


1 


- 5 


0.3 



Table 8: Values of mAOt from a 16^ x 32 lattice at of 5.70 using 219 configurations. 



34 



k 


sink 




t 




0.1600 





0.6887 ±0.0011 


11 - 


14 


0.0 




1 


0.6884 ±0.0012 


11 - 


14 


0.1 




2 


0.6894 ±0.0016 


9 - 


12 


0.2 




3 


0.6900 ±0.0023 


9 - 


12 


0.2 




4 


0.6905 ±0.0031 


9 - 


12 


0.1 


0.1650 





0.4576 ±0.0013 


10 - 


14 


0.2 




1 


0.4574 ±0.0012 


10 - 


13 


0.0 




2 


0.4577 ±0.0014 


9 - 


12 


0.0 




3 


0.4582 ±0.0019 


9 - 


12 


0.1 




4 


0.4589 ±0.0022 


9 - 


12 


0.1 


0.1663 





0.3799 ±0.0016 


8 - 


11 


0.0 




1 


0.3800 ±0.0016 


8 - 


11 


0.1 




2 


0.3822 ±0.0019 


11 - 


14 


0.1 




3 


0.3825 ± 0.0022 


11 - 


14 


0.2 




4 


0.3829 ±0.0026 


9 - 


14 


0.3 


0.1675 





0.2943 ±0.0021 


9 - 


12 


0.1 




1 


0.2946 ±0.0023 


9 - 


12 


0.0 




2 


0.2949 ±0.0023 


9 - 


12 


0.1 




3 


0.2953 ±0.0019 


10 - 


13 


0.1 




4 


0.2955 ± 0.0024 


10 - 


13 


0.1 



Table 9: Values of rriT^a from a 24^ x 32 lattice at jS of 5.70 using 92 configurations. 



35 



k 


sink 






t 


A 


0.1600 





8062 ± 0018 


9 ■ 


- 14 


0.3 




1 


8046 ± 0023 


11 


- 14 


0.2 




2 


8058 ± 0028 


10 


- 13 


0.7 




3 


8076 ± 0047 


9 ■ 


- 12 


0.8 




4 


8022 ± 0056 


6 ■ 


- 10 


0.7 


0.1650 


012 


6601 ± 0031 


6 ■ 


- 11 


1.3 







0.6585 ± 0.0036 


9 ■ 


- 12 


0.3 




1 


6596 ± 0045 


9 ■ 


- 12 


0.5 




2 


6514 ± 0050 


7 . 


- 10 


1.5 




3 


6511 ± 0071 


6 


- 9 


0.8 




4 


6491 ± 0073 


6 


- 9 


0.2 


1 668 


012 


n 6225 + n 0042 


6 


- 9 


6 







0.6246 ±0.0051 


6 


- 9 


0.8 




1 


0.6213 ±0.0064 


7 ■ 


- 10 


0.2 




2 


0.6191 ±0.0084 


7 ■ 


- 10 


0.0 




3 


0.6191 ±0.0086 


7 ■ 


- 10 


0.0 




4 


0.6206 ±0.0103 


7 ■ 


- 10 


0.0 


0.1675 


012 


0.5867 ±0.0062 


6 


- 9 


0.3 







0.5875 ±0.0068 


6 


- 9 


0.7 




1 


0.5863 ±0.0089 


6 


- 9 


0.4 




2 


0.5842 ±0.0117 


6 


- 9 


0.2 




3 


0.5822 ±0.0110 


6 


- 9 


0.3 




4 


0.5912 ±0.0125 


4 


- 8 


0.2 



Table 10: Values of rUpa from a 24^ x 32 lattice at /? of 5.70 using 92 configurations. 



36 



k 


sink 






t 




0.1600 





1.3068 ±0.0051 


8 ■ 


- 11 


0.5 




1 


1.3077 ±0.0063 


9 ■ 


- 12 


0.6 




2 


1.3062 ±0.0065 


6 


- 9 


0.2 




3 


1.3088 ±0.0091 


7 ■ 


- 10 


0.3 




4 


1.3124 ±0.0135 


7- 


- 10 


0.4 


0.1650 


012 


1.0337 ±0.0067 


7- 


- 11 


0.5 







1.0358 ±0.0083 


6 


- 9 


0.1 




1 


1.0313 ±0.0098 


8 ■ 


- 11 


0.2 




2 


1.0295 ±0.0102 


8 ■ 


- 11 


0.1 




3 


1.0323 ±0.0098 


4 


- 7 


0.1 




4 


1.0301 ±0.0104 


4 


- 7 


0.2 


0.1663 


012 


0.9402 ±0.0117 


7- 


- 10 


0.6 







0.9488 ±0.0109 


6 


- 9 


0.4 




1 


0.9440 ±0.0114 


6 


- 9 


0.0 




2 


0.9459 ±0.0118 


5 


- 8 


0.1 




3 


0.9472 ± 0.0135 


4 


- 8 


0.1 




4 


0.9421 ±0.0131 


3 


- 8 


0.1 


0.1675 


012 


0.8634 ±0.0135 


6 


- 9 


1.3 







0.8571 ±0.0118 


6 


- 9 


0.0 




1 


0.8527 ±0.0137 


6 


- 9 


0.2 




2 


0.8579 ±0.0133 


5 


- 8 


0.0 




3 


0.8587 ±0.0158 


5 


- 8 


0.1 




4 


0.8668 ±0.0177 


4 


- 8 


0.1 



Table 11: Values of rriNa from a 24^ x 32 lattice at /? of 5.70 using 92 configurations. 



37 







777 A n 




u 


A 


0.1600 





1 3799 ± 0062 


8 ■ 


- 11 


0.2 




1 


1 Qci 5 + 0061 


9 ■ 


- 12 


2 




2 


1 3826 + 0066 




- 8 









1 3827 + 01 04 




- 8 


1 




4 


1 QQO^ + n 01 63 


7 . 


- 10 


3 


1 fi'iO 


019 




6 ■ 


- 10 


n 6 







1 1614 ± 0101 


7 . 


- 10 


0.6 




1 


1 1693 ± 0128 


5 


- 8 


0.6 




2 


1 1 673 + 0067 


4 


_ 7 


1 

\J . X 






1 1633 + 0107 




- 8 







4 


1 1524 ± 0150 


4 


. 7 


0.1 






1 1 1 70 + n ni 38 


6 . 

U 


- 1 n 


n 7 







1.1120 ± 0.0163 


7 ■ 


- 10 


0.4 




1 


1.1226 ±0.0119 


5 


- 8 


0.7 




2 


1.1170 ±0.0115 


4 


- 8 


0.8 




3 


1.1032 ±0.0164 


3 


- 6 


0.1 




4 


1.0842 ±0.0177 


3 


- 6 


0.0 


0.1675 


012 


1.0817 ±0.0213 


6 


- 9 


1.0 







1.0917± 0.0132 


5 


- 9 


1.1 




1 


1.0783 ±0.0178 


5 


- 9 


0.9 




2 


1.0770 ±0.0128 


4 


- 7 


0.5 




3 


1.0695 ±0.0198 


3 


- 7 


0.1 




4 


1.0523 ±0.0214 


4 


- 7 


0.0 



Table 12: Values of mAOt from a 24^ x 32 lattice at P of 5.70 using 92 configurations. 



38 





olllJv 




t 




X 


n 1 'i4'^ 

yj . ±0^0 


n 




19 - 


1 6 


n 8 

u. 




1 
± 


zn u.uuzu 


j.^ - 


1 7 


1 n 

l.U 




9 
z 


u.^oou zn u.uuzi 


1 1 


1 4 


n 4 




(J 




11 - 


14 


7 

yj. 1 








1 1 - 


14 


1 n 


u.ioou 


n 

U 


u.oo ( u zn u.uuio 


1 9 


1 ^ 
10 


n 4 




1 

X 


fi Qc;cQ -1- n nni 7 

u.oooo zn u.uui / 


1 1 


1 4 


n Q 




9 


fi Qcisn -1- n nni r 

u.ouou zn u.uu±L) 


1 1 


1 4 


n 4 




(J 


3576 zt 0017 


11 - 


14 


3 




4 


0.3573 ±0.0019 


11 - 


14 


0.4 


0.1573 





0.2653 ±0.0017 


12 - 


15 


0.0 




1 


0.2658 ±0.0019 


10 - 


14 


0.2 




2 


0.2653 ±0.0018 


10 - 


13 


0.1 




3 


0.2648 ±0.0020 


10 - 


13 


0.0 




4 


0.2641 ±0.0025 


10 - 


13 


0.0 


0.1581 





0.1906 ±0.0023 


12 - 


15 


0.1 




1 


0.1902 ±0.0023 


11 - 


14 


0.0 




2 


0.1898 ±0.0024 


10 - 


13 


0.0 




3 


0.1892 ±0.0026 


10 - 


13 


0.0 




4 


0.1885 ±0.0031 


10 - 


13 


0.0 



Table 13: Values of m-,^a from a 24^ x 36 lattice at /3 of 5.93 using 210 configurations. 



39 



k 


sink 




t 


x' 


0.1543 





0.5505 ±0.0029 


12 - 


16 


0.8 




1 


0.5528 ±0.0024 


11 - 


14 


0.5 




2 


0.5519 ±0.0027 


11 - 


14 


0.4 




3 


0.5565 ±0.0041 


8 - 


11 


0.4 




4 


0.5527 ±0.0040 


7- 


10 


0.0 


0.1560 


012 


0.4843 ± 0.0033 


12 - 


15 


1.0 







0.4829 ± 0.0042 


12 - 


15 


1.3 




1 


0.4889 ±0.0033 


9 - 


12 


0.9 




2 


0.4846 ± 0.0036 


11 - 


14 


0.3 




3 


0.4891 ± 0.0045 


8 - 


11 


0.1 




4 


0.4864 ± 0.0042 


7- 


10 


0.0 


0.1573 


012 


0.4399 ±0.0047 


9 - 


12 


0.3 







0.4397 ±0.0044 


9 - 


12 


0.6 




1 


0.4393 ± 0.0047 


9 - 


12 


0.5 




2 


0.4404 ± 0.0048 


8 - 


11 


0.3 




3 


0.4398 ± 0.0045 


7- 


10 


0.0 




4 


0.4369 ± 0.0048 


6 - 


9 


0.0 


0.1581 


012 


0.4091 ± 0.0063 


9 - 


12 


0.5 







0.4110 ±0.0058 


9 - 


12 


0.7 




1 


0.4154 ±0.0057 


8 - 


11 


0.6 




2 


0.4146 ±0.0045 


7- 


10 


0.0 




3 


0.4121 ± 0.0042 


7- 


10 


0.2 




4 


0.4071 ± 0.0057 


6 - 


9 


0.1 



Table 14: Values of nipa from a 24^ x 36 lattice at (3 of 5.93 using 210 configurations. 



40 



k 


sink 


rriNCi 




t 




0.1543 





0.8628 ±0.0056 


13 


- 16 


0.6 




1 


0.8662 ±0.0054 


11 


- 14 


0.0 




2 


0.8631 ±0.0048 


11 


- 14 


0.9 




3 


0.8792 ±0.0100 


8 


- 11 


0.6 




4 


0.8674 ±0.0102 


6 


- 9 


0.1 


0.1560 


012 


0.7438 ±0.0053 


12 


- 15 


1.0 







0.7434 ±0.0064 


11 


- 14 


0.4 




1 


0.7427 ±0.0059 


11 


- 14 


0.0 




2 


0.7404 ±0.0062 


11 


- 14 


0.9 




3 


0.7514 ±0.0064 


7 ■ 


- 10 


0.1 




4 


0.7448 ±0.0099 


1 ■ 


- 10 


0.1 


0.1573 


012 


0.6398 ±0.0068 


11 


- 14 


0.9 







0.6391 ± 0.0087 


11 


- 14 


0.3 




1 


0.6389 ±0.0082 


11 


- 14 


0.2 




2 


0.6489 ± 0.0079 


8 ■ 


- 11 


0.1 




3 


0.6469 ± 0.0063 


1 ■ 


- 10 


0.2 




4 


0.6423 ±0.0080 


6 


- 10 


0.1 


0.1581 


012 


0.5759 ±0.0099 


9 


- 12 


1.8 







0.5718 ±0.0074 


9 


- 12 


0.3 




1 


0.5707 ±0.0077 


9 


- 12 


0.1 




2 


0.5714 ±0.0068 


8 ■ 


- 11 


0.1 




3 


0.5701 ± 0.0086 


1 ■ 


- 10 


0.5 




4 


0.5652 ± 0.0092 


6 


- 10 


0.2 



Table 15: Values oirrij^a from a 24^ x 36 lattice at /5 of 5.93 using 210 configurations. 



41 



k 


sink 






t 


x' 


0.1543 





0.9216 ±0.0071 


13 


- 16 


0.3 




1 


0.9280 ±0.0080 


11 


- 14 


0.4 




2 


0.9240 ±0.0076 


11 


- 14 


0.7 




3 


0.9406 ±0.0088 


8 ■ 


- 11 


0.4 




4 


0.9276 ±0.0117 


6 


- 9 


0.1 


0.1560 


012 


0.8314 ±0.0062 


12 


- 15 


1.1 







0.8304 ±0.0085 


11 


- 15 


1.1 




1 


0.8264 ±0.0082 


11 


- 15 


0.5 




2 


0.8225 ±0.0082 


11 


- 14 


0.3 




3 


0.8349 ±0.0081 


7- 


- 10 


0.0 




4 


0.8272 ±0.0115 


7- 


- 10 


0.1 


0.1573 


012 


0.7564 ±0.0103 


11 


- 14 


1.5 







0.7594 ± 0.0098 


10 


- 14 


1.0 




1 


0.7566 ±0.0106 


10 


- 13 


0.5 




2 


0.7620 ±0.0100 


8 ■ 


- 11 


0.2 




3 


0.7581 ± 0.0097 


7- 


- 10 


0.0 




4 


0.7500 ±0.0108 


7- 


- 10 


0.0 


0.1581 


012 


0.7234 ±0.0114 


9 ■ 


- 12 


1.2 







0.7161 ± 0.0094 


9 ■ 


- 12 


0.3 




1 


0.7126 ±0.0113 


9 ■ 


- 12 


0.4 




2 


0.7105 ±0.0113 


8 ■ 


- 11 


0.2 




3 


0.7083 ±0.0114 


7- 


- 10 


0.2 




4 


0.6978 ±0.0141 


7- 


- 10 


0.1 



Table 16: Values of from a 24^ x 36 lattice at /? of 5.93 using 210 configurations. 
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If 


Gm It 




t 




A 


1 ^nn 


n 




14 - 


1 7 


n 6 




1 

X 


U.OO/ u zc u.uuxu 


14 - 


1 7 


1 n 




9 


n '^SR'^ -1- OOOQ 


±"4: - 


1 Q 


1 9 




3 




13 - 


18 


Q 




4 




13 - 


1 8 


n 7 


n 1 ^1 Q 


n 


n ORQQ -L n nni o 


14 - 


1 7 


n 3 




1 

± 


n ORQQ + n nni o 


14 - 


1 7 


n 7 

yj. 1 




9 


n 9«Qn -L n nni o 

U.ZiUOU ZIZ U.UUXZi 


14 - 




R 
u.u 




3 


2631 ± 0011 


14 - 


18 


0.4 




4 


0.2631 ±0.0012 


14 - 


18 


0.3 


0.1526 





0.2065 ±0.0017 


14 - 


17 


0.4 




1 


0.2066 ±0.0016 


14 - 


17 


0.3 




2 


0.2066 ±0.0014 


14 - 


18 


0.2 




3 


0.2067 ±0.0014 


14 - 


18 


0.1 




4 


0.2064 ±0.0015 


15 - 


18 


0.1 


0.1532 





0.1445 ±0.0020 


12 - 


15 


0.2 




1 


0.1453 ±0.0020 


15 - 


18 


0.4 




2 


0.1451 ±0.0019 


13 - 


17 


0.1 




3 


0.1453 ±0.0020 


13 - 


16 


0.0 




4 


0.1455 ±0.0020 


13 - 


16 


0.0 



Table 17: Values of m^a from a 32 x 30 x 32 x 40 lattice at (3 of 6.17 using 219 
configurations. 
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k 


sink 




t 






0.1500 





0.4462 ±0.0012 


14 - 


17 


0.1 




1 


0.4461 ±0.0011 


14 - 


17 


0.2 




2 


0.4458 ±0.0012 


12 - 


16 


0.6 




3 


0.4458 ±0.0014 


13 - 


16 


0.5 




4 


0.4458 ±0.0018 


13 - 


16 


0.4 


0.1519 


012 


0.3581 ±0.0018 


13 - 


16 


0.3 







0.3577 ±0.0017 


14 - 


17 


0.1 




1 


0.3577 ±0.0018 


13 - 


16 


0.2 




2 


0.3573 ±0.0018 


12 - 


16 


0.2 




3 


0.3573 ±0.0022 


13 - 


16 


0.3 




4 


0.3572 ± 0.0026 


13 - 


16 


0.3 


0.1526 


012 


0.3265 ± 0.0022 


13 - 


16 


0.4 







0.3254 ±0.0027 


13 - 


16 


0.1 




1 


0.3252 ±0.0026 


13 - 


16 


0.0 




2 


0.3247 ±0.0027 


12 - 


16 


0.1 




3 


0.3246 ±0.0032 


13 - 


16 


0.1 




4 


0.3245 ± 0.0039 


13 - 


16 


0.1 


0.1532 


012 


0.3024 ± 0.0048 


13 - 


16 


0.4 







0.2991 ± 0.0059 


13 - 


16 


0.2 




1 


0.2987 ±0.0059 


12 - 


15 


0.1 




2 


0.2969 ±0.0066 


13 - 


16 


0.1 




3 


0.2965 ± 0.0076 


13 - 


16 


0.0 




4 


0.2965 ± 0.0088 


13 - 


16 


0.0 



Table 18: Values of m^a from a 32 x 30 x 32 x 40 lattice at (3 of 6.17 using 219 
configurations. 
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k 


sink 




t 




x' 


0.1500 





0.7004 ±0.0028 


14 - 


17 


0.3 




1 


0.6995 ±0.0025 


14 - 


17 


0.7 




2 


0.6975 ±0.0029 


14 - 


18 


0.9 




3 


0.6970 ±0.0030 


14 - 


18 


0.7 




4 


0.6966 ±0.0040 


14 - 


18 


0.6 


0.1519 


012 


0.5427 ±0.0038 


15 - 


18 


0.9 







0.5444 ±0.0037 


14 - 


17 


0.2 




1 


0.5437 ±0.0035 


14 - 


17 


0.8 




2 


0.5424 ±0.0036 


14 - 


18 


1.0 




3 


0.5428 ±0.0038 


14 - 


18 


0.6 




4 


0.5460 ±0.0052 


13 - 


17 


0.4 


0.1526 


012 


0.4845 ± 0.0053 


13 - 


16 


1.5 







0.4819 ±0.0050 


14 - 


17 


0.9 




1 


0.4826 ± 0.0048 


12 - 


16 


0.6 




2 


0.4824 ±0.0050 


13 - 


16 


0.3 




3 


0.4831 ± 0.0042 


13 - 


16 


0.2 




4 


0.4848 ± 0.0068 


13 - 


17 


0.1 


0.1532 


012 


0.4234 ± 0.0065 


13 - 


16 


1.1 







0.4187 ±0.0059 


12 - 


15 


0.5 




1 


0.4188 ±0.0087 


12 - 


17 


0.7 




2 


0.4199 ±0.0088 


13 - 


16 


0.5 




3 


0.4215 ±0.0080 


13 - 


16 


0.2 




4 


0.4097 ±0.0078 


10 - 


13 


0.9 



Table 19: Values of mjva from a 32 x 30 x 32 x 40 lattice at (3 of 6.17 using 219 
configurations. 
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k 


sink 




t 




x' 


0.1500 





0.7406 ±0.0034 


14 - 


17 


0.5 




1 


0.7398 ±0.0030 


14 - 


17 


0.7 




2 


0.7379 ±0.0036 


14 - 


18 


0.6 




3 


0.7368 ±0.0038 


14 - 


18 


0.4 




4 


0.7358 ±0.0045 


14 - 


18 


0.3 


0.1519 


012 


0.6126 ±0.0041 


13 - 


17 


0.6 







0.6107 ±0.0040 


14 - 


17 


0.2 




1 


0.6101 ±0.0039 


14 - 


17 


0.7 




2 


0.6087 ±0.0047 


12 - 


16 


0.5 




3 


0.6085 ± 0.0046 


13 - 


16 


0.5 




4 


0.6083 ±0.0058 


13 - 


16 


0.4 


0.1526 


012 


0.5656 ±0.0046 


13 - 


17 


0.4 







0.5637 ±0.0049 


14 - 


17 


0.0 




1 


0.5626 ±0.0053 


14 - 


17 


0.2 




2 


0.5610 ±0.0056 


12 - 


16 


0.5 




3 


0.5608 ±0.0063 


13 - 


16 


0.4 




4 


0.5598 ±0.0073 


13 - 


16 


0.4 


0.1532 


012 


0.5267 ±0.0081 


13 - 


17 


0.9 







0.5233 ±0.0088 


14 - 


17 


0.1 




1 


0.5188 ±0.0111 


14 - 


17 


0.0 




2 


0.5133 ±0.0104 


14 - 


17 


0.5 




3 


0.5197 ±0.0101 


11 - 


16 


0.9 




4 


0.5120 ±0.0101 


11 - 


15 


0.7 



Table 20: Values of niAa from a 32 x 30 x 32 x 40 lattice at P of 6.17 using 219 
configurations. 
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lattice 


/3 


sink 






rUpa 


16=^ X 32 


5.70 


012 


0.00390 ± 0.00012 


0.09662 ± 0.00291 


0.5676 ± 0.0079 






4 


0.00372 ± 0.00018 


0.09223 ± 0.00434 


0.5522 ±0.0132 


24^ X 32 


5.70 


012 


0.00348 ± 0.00012 


0.08620 ± 0.00303 


0.5409 ±0.0089 






4 


0.00366 ±0.00029 


0.09070 ±0.00710 


0.5560 ±0.0212 


243 X 36 


5.93 


012 


0.00223 ±0.00009 


0.05536 ±0.00226 


0.3851 ±0.0079 






4 


0.00216 ±0.00008 


0.05364 ±0.00191 


0.3803 ±0.0064 


30 X 322 X 4Q 


6.17 


012 


0.00141 ±0.00004 


0.03503 ±0.00096 


0.2768 ±0.0039 






4 


0.00139 ± 0.00007 


0.03451 ± 0.00164 


0.2736 ± 0.0065 



Table 21: The normal quark mass m^a and strange quark mass mgO, and the rho 
mass extrapolated to quark mass m„a. 
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sink 






rriK* 1 TUp 


012 


1.106 ±0.009 


0.6 




4 


1.123 d= 0.015 


0.7 




012 


1.213 ±0.017 


0.6 




4 


1.246 ±0.030 


0.7 


rriN/mp 


012 


1.459 ±0.034 


1.3 




4 


1.371 ±0.038 


1.9 


/\m/mp 


012 


1.914 ±0.031 


1.3 




4 


1.939 ±0.027 


1.9 


m 1^1 nip 


012 


1.922 ±0.036 


3.4 




4 


1.786 ± 0.048 


0.0 


niY,*/nip 


012 


1.986 ±0.024 


3.4 




4 


1.903 ±0.032 


0.0 


TTis* 1 Tap 


012 


2.049 ±0.030 


3.4 




4 


2.020 ±0.037 


0.0 


m^jnip 


012 


2.112 ±0.050 


3.4 




4 


2.137 ±0.043 


0.0 




012 


0.301 ± 0.007 




4 


0.301 ± 0.007 





Table 22: Hadron masses extrapolated to physical quark mass measured in units of 
nip extrapolated to physical quark mass from a 16^ x 32 lattice at (3 of 5.70. The 
mass difference Am is ms ± m^, — m^- 
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sink 






rriK* 1 Tap 


012 


1.124 ±0.007 


0.0 




4 


1.099 ±0.023 


0.0 




012 


1.249 ±0.014 


0.0 




4 


1.198 ±0.046 


0.0 


rriN/mp 


012 


1.397 ±0.036 


0.2 




4 


1.373 ±0.064 


0.0 


/\m/mp 


012 


1.975 ±0.015 


0.2 




4 


1.940 ±0.035 


0.0 


m 1^1 nip 


012 


1.879 ±0.052 


0.4 




4 


1.775 ± 0.090 


0.9 


niY,*/nip 


012 


1.987 ±0.033 


0.4 




4 


1.890 ± 0.070 


0.9 


TTis* 1 Tap 


012 


2.095 ±0.022 


0.4 




4 


2.005 ± 0.056 


0.9 


m^jnip 


012 


2.203 ±0.028 


0.4 




4 


2.119 ±0.055 


0.9 




012 


0.299 ±0.011 




4 


0.299 ±0.011 





Table 23: Hadron masses extrapolated to physical quark mass measured in units of 
nip extrapolated to physical quark mass from a 24^ x 32 lattice at (3 of 5.70. The 
mass difference Am is ms ± m^, — m^- 
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sink 












rriK* 1 TUp 


012 


1, 


.121 ± 


0, 


.008 


0.4 




4 


1, 


.126 ± 


0, 


.007 


0.0 




012 


1, 


.242 ± 


0, 


.016 


0.4 




4 


1, 


.252 ± 


0, 


.013 


0.0 


rriN/mp 


012 


1, 


.347 ± 


0, 


.028 


0.0 




4 


1 


.340 ± 


0, 


.031 


1.9 


/\m/mp 


012 


1, 


.890 ± 


0, 


.014 


0.0 




4 


1, 


.905 ± 


0, 


.027 


1.9 


m 1^1 nip 


012 


1 


.772 ± 


0, 


.038 


0.6 




4 


1 


.728 ± 


0, 


.036 


0.2 


niY,*/nip 


012 


1 


.893 ± 


0, 


.034 


0.6 




4 


1 


.863 ± 


0, 


.033 


0.2 


TTis* 1 Tap 


012 


2, 


.014 ± 


0, 


.030 


0.6 




4 


1 


.998 ± 


0, 


.029 


0.2 


m^jnip 


012 


2, 


.135 ± 


0, 


.031 


0.6 




4 


2, 


.134 ± 


0, 


.035 


0.2 




012 


0, 


.301 ± 


0, 


.005 




4 


0, 


.301 ± 


0, 


.005 





Table 24: Hadron masses extrapolated to physical quark mass measured in units of 
nip extrapolated to physical quark mass from a 24^ x 36 lattice at (3 of 5.93. The 
mass difference Am is ms + m^, — m^- 
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sink 












rriK* 1 TUp 


012 


1, 


.128 ± 


0, 


.004 


0.6 




4 


1, 


.131 ± 





.004 


0.0 




012 


1, 


.256 ± 


0, 


.007 


0.6 




4 


1, 


.262 ± 


0, 


.009 


0.0 


rriN/mp 


012 


1, 


.383 ± 


0, 


.030 


2.2 




4 


1 


.360 ± 


0, 


.040 


4.6 


/\m/mp 


012 


1, 


.888 ± 


0, 


.015 


2.2 




4 


1, 


.886 ± 


0, 


.024 


4.6 


m 1^1 nip 


012 


1 


.782 ± 


0, 


.039 


0.0 




4 


1 


.751 ± 


0, 


.055 


0.3 


niY,*/nip 


012 


1 


.907 ± 


0, 


.031 


0.0 




4 


1 


.885 ± 


0, 


.046 


0.3 


TTis* 1 Tap 


012 


2, 


.033 ± 


0, 


.025 


0.0 




4 


2, 


.018 ± 


0, 


.037 


0.3 


m^jnip 


012 


2, 


.158 ± 


0, 


.018 


0.0 




4 


2, 


.152 ± 


0, 


.032 


0.3 




012 


0, 


.302 ± 


0, 


.007 




4 


0, 


.302 ± 


0, 


.007 





Table 25: Hadron masses extrapolated to physical quark mass measured in units of 
nip extrapolated to physical quark mass from a 30 x 32^ x 40 lattice at /3 of 6.17. The 
mass difference Am is ms + ms — m^- 
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sink 






fin. vol. 




inf. vol. 




obs. 


rriK* 1 mp 


012 


0.01 


1, 


.149±0.009 


1, 


.166±0.016 


1, 


.164 




4 


0.03 


1, 


.140±0.016 


1, 


.116±0.032 








012 


0.01 


1, 


.297±0.019 


1, 


.333±0.032 


1, 


.327 




4 


0.03 


1, 


.280±0.032 


1, 


.232±0.064 






rriN/mp 


012 


2.61 


1, 


.278±0.068 


1, 


.216±0.104 


1, 


.222 




4 


0.27 


1, 


.328±0.085 


1, 


.330±0.131 






Am /nip 


012 


0.12 


L 


.865±0.045 


1, 


.927±0.074 


2, 


.047 




4 


0.00 


1, 


.834±0.055 


1, 


.835±0.083 






rriA/mp 


012 


1.37 


1 


.609±0.078 


1 


.565±0.122 


1 


.604 




4 


0.42 


1 


.683±0.094 


1 


.672±0.174 






mY,*/mp 


012 


1.04 


1 


.805±0.058 


1 


.806±0.080 


1 


.803 




4 


0.37 


1 


.840±0.084 


1 


.827±0.136 






rriB,* 1 mp 


012 


0.48 


2, 


.009±0.049 


2, 


.055±0.065 


1 


.996 




4 


0.23 


2, 


.003±0.079 


1 


.988±0.111 






m^/mp 


012 


0.03 


2, 


.205±0.056 


2, 


.296±0.089 


2, 


.177 




4 


0.08 


2, 


.164±0.072 


2, 


.147±0.121 








012 


0.02 


0, 


.305±0.008 


0, 


.319±0.012 


0, 


.305 ± 0.018 


4 


0.01 


0, 


.303±0.012 


0, 


.301±0.020 







Table 26: Calculated values of hadron mass ratios at physical quark masses, ex- 
trapolated to zero lattice spacing in finite volume, then corrected to infinite volume, 
compared with observed values. The mass difference Am is tos + ms — m^- 
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k 


sink 012 


sink 4 


rrin 


U.iooU 
0.1663 
0.1675 


— U.2 ± U.o/o 
-1.8 ±0.7% 
-1.1 ± 1.5% 


u.2 ± U.o/o 
-0.9 ± 1.1% 
-1.2 ± 1.9% 


rup 


U.iDOU 

0.1663 
0.1675 


— U.D ± u.y/o 
-1.5 ± 1.4% 
-3.4 ± 1.4% 


— o.o =t i.o/o 
-0.9 ±2.1% 
-1.6 ±2.6% 


rriN 


0.1650 
0.1663 
0.1675 


0.2 ± 1.4% 
-4.6 ± 1.8% 
-4.8 ± 2.0% 


-0.2 ± 1.8% 
-1.6 ±2.0% 
0.7 ± 2.4% 


ruA 


0.1650 
0.1663 

0.1675 


0.1 ± 1.6% 
-0.4 ± 2.0% 

-1.9± 2.2'X 


0.1 ±2.1% 
-1.6 ±2.6% 

-0.0 ± 2.3% 



Table 27: Changes in masses from 16^ x 32 to 24^ x 32 at /? = 5.70. 





sink 012 


sink 4 




1.6±U.9% 


-2.2 ± 2.2% 


m<p/mp 


2.9 ± 1.6% 


-4.0 ± 4.U% 


rriN/mp 


-4.5 ± 3.7% 


U.2 ± 5.8% 


t^mjmp 


3.1 ± 2.1% 


0.0 ±2.5% 


rriA/mp 


-2.3 ±3.3% 


-0.6 ± 5.8% 


mj:*/mp 


0.1 ±2.1% 


-0.7 ± 4.4% 


ms*/mp 


2.2 ± 1.7% 


-0.8 ± 3.3% 


ma/ nip 


4.1 ± 2.4% 


-0.8 ± 3.2% 




4.7 ±2.1% 


-0.7 ±4.3% 



Table 28: Changes in mass ratios from 16^ x 32 to 24^ x 32 at /3 = 5.70 The mass 
difference Am is ms + ms — tun- 
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sink 








mp{m„s) druq ^ 


012 


0.093 ±0, 


.008 


0.6 




4 


0.109 ±0, 


.015 


0.7 


mjsi{mns) 


012 


1.517±0, 


.014 


1.3 




4 


1.465 ±0, 


.022 


1.9 


mp{inns) drriq ^ 


012 


0.199 ±0, 


.026 


1.3 




4 


0.253 ±0, 


.017 


1.9 


mA(m„s) 
mp{mns) 


012 


1.834 dzO, 


.016 


3.4 




4 


1.750 ±0, 


.024 


0.0 


rUna dm A /jj^ \ 
mp{m„s) dmq V 


012 


0.083 ±0, 


.029 


3.4 




4 


0.156 ±0, 


.021 


0.0 


a(0) 










MS 

mp{mns) 


012 


0.2677 ± 0. 


,0027 






4 


0.2745 ± 0. 


,0039 





Table 29: Hadron mass parameters determined without extrapolation in quark mass 
for the lattice 16^ x 32 lattice at (5 of 5.70. 
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sink 








mp{m„s) drriq ^ 


012 


0.111 ±0, 


.007 


0.0 




4 


0.086 ±0, 


.022 


0.0 


m]si{mns) 


012 


1.492 ±0, 


.017 


0.2 




4 


1.490 ±0, 


.026 


0.0 


mp{m„s) drriq ^ 


012 


0.258 ±0, 


.017 


0.2 




4 


0.246 ±0, 


.027 


0.0 


mA(m„s) 
mp{mns) 


012 


1.820 ±0, 


.023 


0.4 




4 


1.772 lb 0, 


.039 


0.9 


m„a dm A /jj^ \ 
mp{m„s) dmq V 


012 


0.144 ±0, 


.027 


0.4 




4 


0.149 ±0, 


.033 


0.9 


a(0) 










MS 

mp{mns) 


012 


0.2759 ± 0. 


,0029 






4 


0.2752 ± 0. 


,0047 





Table 30: Hadron mass parameters determined without extrapolation in quark mass 
for the lattice 24^ x 32 lattice at (5 of 5.70. 
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sink 








mp{m„s) drriq ^ 


012 


0.107 ±0, 


.008 


0.4 




4 


0.113 ±0, 


.007 


0.0 


m]si{mns) 


012 


1.435 ±0, 


.015 


0.0 




4 


1.433 ±0, 


.017 


1.9 


mp{m„s) drriq ^ 


012 


0.241 ±0, 


.010 


0.0 




4 


0.252 ±0, 


.018 


1.9 


mA(m„s) 
mp{mns) 


012 


1.745 ±0, 


.023 


0.6 




4 


1.715 ±0, 


.024 


0.2 


m„a dm A /jj^ \ 
mp{m„s) dmq V 


012 


0.161 ±0, 


.012 


0.6 




4 


0.181 ±0, 


.012 


0.2 


a(0) 










MS 

mp{mns) 


012 


0.2679 ± 0. 


,0041 






4 


0.2698 ± 0. 


,0034 





Table 31: Hadron mass parameters determined without extrapolation in quark mass 
for the lattice 24^ x 36 lattice at (5 of 5.93. 
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sink 






mp{m„s) druq ^ 


012 


0.115 ±0.004 


0.6 




4 


0.118 ±0.005 


0.0 


mjsi{mns) 


012 


1.445 ±0.017 


2.2 




4 


1.430 ±0.024 


4.6 


mp{inns) drriq ^ 


012 


0.226 ±0.010 


2.2 




4 


0.236 ±0.012 


4.6 


mA(m„s) 
mp{mns) 


012 


1.749 ±0.019 


0.0 




4 


1.727 ±0.033 


0.3 


rUna dm A /jj^ \ 
mp{m„s) dmq V 


012 


0.169 ±0.011 


0.0 




4 


0.180 ±0.012 


0.3 


a(0) 








MS 

mp{mns) 


012 


0.2671 ± 0.0032 






4 


0.2695 ± 0.0057 





Table 32: Hadron mass parameters determined without extrapolation in quark mass 
for the lattice 30 x 32^ x 40 lattice at /3 of 6.17. 
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sink 




fin. vol. 


inf. vol. 


obs. 


mp{mns) dniq ^ 


012 


0.0 


0.136 ±0.010 


0.154 ±0.016 


0.152 




4 


0.0 


0.127 ±0.017 


0.103 ±0.032 




m,jv(mns) 


012 


3.2 


1.342 ± 0.031 


1.317 ±0.045 


1.405 




4 


0.2 


1.384 ± 0.047 


1.409 ±0.071 




mp{mns) druq ^ 


012 


2.5 


0.239 ±0.028 


0.298 ±0.055 


0.384 




4 


0.2 


0.221 ± 0.028 


0.215 ±0.062 




rnAirrina) 
mp{rans) 


012 


1.6 


1.639 ± 0.040 


1.625 ± 0.053 


1.628 




4 


0.4 


1.684 ±0.053 


1.707 ±0.082 




mp{mna) dniq ^ '"■^) 


012 


1.5 


0.245 ± 0.032 


0.306 ± 0.059 


0.267 


MS 

mp{mns) 


4 


0.2 


0.205 ± 0.041 


0.198 ±0.075 




012 




0.2669 ±0.0057 


0.2750 ±0.0077 


0.262 ±0.015 




4 




0.2698 ± 0.0034 


0.2728 ±0.0120 





Table 33: Calculated values of hadron mass parameters determined without ex- 
trapolation in quark mass, extrapolated to zero lattice spacing in finite volume, then 
corrected to infinite volume, compared with observed values. 
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Figure 1: Fractional uncertainty in the pion propagator obtained from one source 
(circles) and from eight sources (triangles) on the lattice 16^ x 32 at /9 = 5.70 and 
k = 0.1650. 
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Figure 2: Fractional uncertainty in the pion propagator obtained from one source 
(circles) and from eight sources (triangles) on the lattice 16^ x 32 at /9 = 5.70 and 
k = 0.1675. 
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Figure 3: Fractional uncertainty in the nucleon propagator obtained from one source 
(circles) and from eight sources (triangles) on the lattice 16^ x 32 at /9 = 5.70 and 
k = 0.1650. 
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Figure 4: Fractional uncertainty in the nucleon propagator obtained from one source 
(circles) and from eight sources (triangles) on the lattice 16^ x 32 at /9 = 5.70 and 
k = 0.1675. 
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Figure 5: Fractional uncertainty in the pion propagator obtained from one source 
(circles) and from eight sources (triangles) on the lattice 24^ x 32 at /9 = 5.70 and 
k = 0.1650. 
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Figure 6: Fractional uncertainty in the pion propagator obtained from one source 
(circles) and from eight sources (triangles) on the lattice 24^ x 32 at /9 = 5.70 and 
k = 0.1675. 
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Figure 7: Fractional uncertainty in the nucleon propagator obtained from one source 
(circles) and from eight sources (triangles) on the lattice 24^ x 32 at /3 = 5.70 and 
k = 0.1650. 
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Figure 8: Fractional uncertainty in the nucleon propagator obtained from one source 
(circles) and from eight sources (triangles) on the lattice 24^ x 32 at /3 = 5.70 and 
k = 0.1675. 
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Figure 9: Effective masses, final fitting range and fitted mass for tfie pseudoscalar 
propagator with sink size 2 on tfie lattice 16^ x 32 at at /3 = 5.70 and k — 0.1675 



67 



0.40 



0.35 



E 0.30 



0.25 



0.20 



t 







8 



12 



16 



t 



Figure 10: Effective masses, final fitting range and fitted mass for tfie pseudoscalar 
propagator with sink size 4 on tfie lattice 16^ x 32 at at /3 = 5.70 and k — 0.1675 
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Figure 11: Effective masses, final fitting range and fitted mass for tfie pseudoscalar 
propagator with sink size 2 on tfie lattice 24^ x 36 at at /3 = 5.93 and k — 0.1581 
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Figure 12: Effective masses, final fitting range and fitted mass for tfie pseudoscalar 
propagator with sink size 4 on tfie lattice 24^ x 36 at at /3 = 5.93 and k — 0.1581 
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Figure 13: Effective masses, final fitting range and fitted mass for tfie pseudoscalar 
propagator with sink size 2 on tfie lattice 30 x 32^ x 40 at at /3 = 6.17 and k — 0.1532 
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Figure 14: Effective masses, final fitting range and fitted mass for tfie pseudoscalar 
propagator with sink size 4 on tfie lattice 30 x 32^ x 40 at at /3 = 6.17 and k — 0.1532 
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Figure 15: Effective masses, final fitting range and fitted mass for the vector propa- 
gator with sink size 2 on the lattice 16^ x 32 at at /3 = 5.70 and k — 0.1675 
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Figure 16: Effective masses, final fitting range and fitted mass for the vector propa- 
gator with sink size 4 on the lattice 16^ x 32 at at /3 = 5.70 and k — 0.1675 
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Figure 17: Effective masses, final fitting range and fitted mass for the vector propa- 
gator with sink size 2 on the lattice 24^ x 36 at at /3 = 5.93 and k — 0.1581 
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Figure 18: Effective masses, final fitting range and fitted mass for the vector propa- 
gator with sink size 4 on the lattice 24^ x 36 at at /3 = 5.93 and k — 0.1581 
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Figure 19: Effective masses, final fitting range and fitted mass for the vector propa- 
gator with sink size 2 on the lattice 30 x 32^ x 40 at at /3 = 6.17 and k = 0.1532 
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Figure 20: Effective masses, final fitting range and fitted mass for the vector propa- 
gator with sink size 4 on the lattice 30 x 32^ x 40 at at /3 = 6.17 and k = 0.1532 
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Figure 21: Effective masses, final fitting range and fitted mass for the nucleon 
propagator with sink size 2 on the lattice 16^ x 32 at at /3 = 5.70 and k = 0.1675 
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Figure 22: Effective masses, final fitting range and fitted mass for the nucleon 
propagator with sink size 4 on the lattice 16^ x 32 at at /3 = 5.70 and k — 0.1675 
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Figure 23: Effective masses, final fitting range and fitted mass for the nucleon 
propagator with sink size 2 on the lattice 24=^ x 36 at at /3 = 5.93 and k = 0.1581 
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Figure 24: Effective masses, final fitting range and fitted mass for the nucleon 
propagator with sink size 4 on the lattice 24=^ x 36 at at /3 = 5.93 and k = 0.1581 
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Figure 25: Effective masses, final fitting range and fitted mass for the nucleon 
propagator with sink size 2 on the lattice 30 x 32^ x 40 at at /3 = 6.17 and k = 0.1532 
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Figure 26: Effective masses, final fitting range and fitted mass for the nucleon 
propagator with sink size 4 on the lattice 30 x 32^ x 40 at at /3 = 6.17 and k = 0.1532 
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Figure 27: Effective masses, final fitting range and fitted mass for the delta baryon 
propagator with sink size 2 on the lattice 16^ x 32 at at /3 = 5.70 and k — 0.1675 
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Figure 28: Effective masses, final fitting range and fitted mass for the delta baryon 
propagator with sink size 4 on the lattice 16^ x 32 at at /3 = 5.70 and k — 0.1675 
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Figure 29: Effective masses, final fitting range and fitted mass for the delta baryon 
propagator with sink size 2 on the lattice 24=^ x 36 at at /3 = 5.93 and k = 0.1581 
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Figure 30: Effective masses, final fitting range and fitted mass for the delta baryon 
propagator with sink size 4 on the lattice 24=^ x 36 at at /3 = 5.93 and k = 0.1581 
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Figure 31: Effective masses, final fitting range and fitted mass for the delta baryon 
propagator with sink size 2 on the lattice 30 x 32^ x 40 at at /3 = 6.17 and k = 0.1532 
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Figure 32: Effective masses, final fitting range and fitted mass for the delta baryon 
propagator with sink size 4 on the lattice 30 x 32^ x 40 at at /3 = 6.17 and k = 0.1532 
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Figure 33: For a 16^ x 32 lattice at (3 of 5.70 combining sinks of sizes 0, 1 and 2, 



m^, nip, rriN and uia, in units of the physical rho mass m. 



as functions of the 



quark mass niq, in units of the strange quark mass m^. The symbol at each point is 
larger than the error bars. 
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Figure 34: For a 16^ x 32 lattice at /3 of 5.70 with a sink of size 4, m^, nip, rriN and 
mA, in units of the physical rho mass mp{mn), as functions of the quark mass rrig, in 
units of the strange quark mass m^. The symbols at most points are larger than the 
error bars. 
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Figure 35: For a 24^ x 36 lattice at (3 of 5.93 combining sinks of sizes 0, 1 and 2, 
m^, nip, rriN and uia, in units of the physical rho mass mp{mn), as functions of the 
quark mass niq, in units of the strange quark mass m^. The symbol at each point is 
larger than the error bars. 
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Figure 36: For a 24^ x 32 lattice at /3 of 5.93 with a sink of size 4, m^, nip, rriN and 
mA, in units of the physical rho mass mp{mn), as functions of the quark mass nig, 
in units of the strange quark mass nis- The symbol at each point is larger than the 
error bars. 
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Figure 37: For a 30 x 32^ x 40 lattice at (3 of 6.17 combining sinks of sizes 0, 1 and 
2, m^, nip, rriN and m^, in units of the physical rho mass mp{mn), as functions of 
the quark mass iriq, in units of the strange quark mass m^. The symbol at each point 
is larger than the error bars. 
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Figure 38: For a 30 x 32^ x 40 lattice at /3 of 6.17 with a sink of size 4, m^, m^, uin 
and TTT-A, in units of the physical rho mass mp{mn)-i as functions of the quark mass 
mq, in units of the strange quark mass m^. The symbol at each point is larger than 
the error bars. 
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Figure 39: Linear extrapolation down to rriq = rUn of synthetic values of mp{mq), 
?7isA(?7T.n) "^q) a^d m/\{rnq) obtained from strange hadron masses. 
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Figure 40: For sinks 0, 1 and 2 combined, mx*/mp and m^jmp as functions of the 
lattice spacing a, in units of l/rrip. The straight hnes are extrapolations to zero lattice 
spacing, the error bars at zero lattice spacing are uncertainties in the extrapolated 
ratios, and the points at zero lattice spacing are observed values. 
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Figure 41: For sink 4, mK*/'^p and m^/mp as functions of the lattice spacing a, 
in units of I /nip. The straight hues are extraploations to zero lattice spacing, the 
error bars at zero lattice spacing are uncertainties in the extrapolated ratios, and the 
points at zero lattice spacing are observed values. 
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Figure 42: For sinks 0, 1 and 2 combined, mjv/mp and (mg + mx; — 111^) /nip as func- 
tions of the lattice spacing a, in units of l/nip. The straight hues are extraploations 
to zero lattice spacing, the error bars at zero lattice spacing are uncertainties in the 
extrapolated ratios, and the points at zero lattice spacing are observed values. 
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Figure 43: For sink 4, m^v/m^ and (m= + — mjv)/'^p as functions of the lattice 
spacing a, in units of l/nip. The straight hues are extraploations to zero lattice 
spacing, the error bars at zero lattice spacing are uncertainties in the extrapolated 
ratios, and the points at zero lattice spacing are observed values. 
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Figure 44: For sinks 0, 1 and 2 combined, m^/nip, 111^,* /irip, m^* /nip and m^/mp as 
functions of the lattice spacing a, in units of l/rrip. The straight hues are extraploat- 
ions to zero lattice spacing, the error bars at zero lattice spacing are uncertainties in 
the extrapolated ratios, and the points at zero lattice spacing are observed values. 
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Figure 45: For sink 4, m^^/mp^ ms*/'mp, m^*/mp and niQ/nip as functions of the 
lattice spacing a, in units of 1/mp. The straight hues are extraploations to zero lattice 
spacing, the error bars at zero lattice spacing are uncertainties in the extrapolated 
ratios, and the points at zero lattice spacing are observed values. 
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Figure 46: For sinks 0, 1 and 2 combined, A^/m^ as a function of the lattice 
spacing a, in units of I /nip. The straight hne is an extraploations to zero lattice 
spacing, the error bars at zero lattice spacing are the uncertainty in the extrapolated 
ratio, and the point at zero lattice spacing is another groups result. 
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Figure 47: For sink 4, A^/rUp as a function of the lattice spacing a, in units of 
l/rUp. The straight hne is an extraploations to zero lattice spacing, the error bars at 
zero lattice spacing are the uncertainty in the extrapolated ratio, and the point at 
zero lattice spacing is another groups result. 
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Figure 48: For sinks 0, 1 and 2, m^a as a function of o;^^ in comparison to the 
prediction of the Callan-Symanzik equation using the two-loop beta function and the 
physical value of found from the continuum hmit of 
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Figure 49: For sinks 0, 1 and 2, hadron masses evaluated at quark mass rrins as a 
function of lattice spacing in units of l/mp{mns)- The straight lines are extraploations 
to zero lattice spacing, the error bars at zero lattice spacing are uncertainties in the 
extrapolated ratios, and the points at zero lattice spacing are observed values. 
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Figure 50: For sinks 0, 1 and 2, hadron masses slopes s defined as dm{mq) / druq 
evaluated at quark mass m„s as a function of lattice spacing in units of l/mp{mns)- 
The straight lines are extraploations to zero lattice spacing, the error bars at zero 
lattice spacing are uncertainties in the extrapolated ratios, and the points at zero 
lattice spacing are observed values. 
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